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ABSTRACT 


An  understanding  of  the  processes  responsible  for  driving  and  damping  acoustic  oscillations 
in  solid  rocket  motors  is  necessary  for  developing  practical  design  methods  that  eliminate  or 
reduce  the  occurrence  combustion  instabilities.  While  state  of  the  art  solid  rocket  stability 
prediction  methods  generally  account  for  the  flow  turning  loss,  the  magnitude  and  characteristics 
of  this  loss  have  never  been  fiilly  investigated.  This  report  describes  results  of  an  investigation  of 
the  role  of  the  flow  turning  loss  in  the  stability  of  solid  rockets  and  its  dependence  upon  motor 
design  and  operating  parameters.  A  one  dimensional  acoustic  stability  equation  that  verifies  that 
the  flow  turning  loss  term  is  appropriately  included  in  the  one  dimensional  stability  formulation 
was  derived  for  a  chamber  with  a  constant  mean  temperature  and  pressure  by  an  approach 
independent  from  that  of  Culick.  This  study  was  extended  provided  the  background  and 
expressions  needed  to  guide  an  experimental  study  of  the  flow  turning  loss  in  the  presence  of 
mean  temperature  and  density  gradients.  This  allows  the  study  of  combustion  systems  in  which 
mean  temperature  gradients  and  heat  losses  are  significant.  The  relevant  conservation  equations 
were  solved  numerically  for  the  experimental  configuration  in  order  to  predict  the  behavior  of  the 
flow  turning  loss  and  to  assist  in  the  analysis  of  experimental  results.  Experiments  performed, 
with  and  without  combustion,  showed  that  the  flow  turning  loss  strongly  depends  upon  the 
propellant  burning  rate  and  the  location  of  the  flow  turning  region  relative  to  the  standing 
pressure  wave.  While  the  stability  equation  for  flows  with  temperature  gradients  differs  fi'om  the 
corresponding  equation  for  cold  flows,  the  measured  behavior  of  the  flow  turning  loss  term  was 
qualitatively  similar  in  experiments  with  and  without  combustion. 
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INTRODUCTION 


The  objective  of  this  research  program  was  to  develop  a  better  understanding  of  the  flow 
turning  loss  and  other  processes  that  contribute  to  the  driving  and  damping  of  axial  instabilities  in 
solid  rocket  motors.  This  understanding  is  necessary  for  the  development  of  dependable  methods 
for  designing  motors  that  do  not  exhibit  such  instabilities. 

The  occurrence  of  combustion  instabilities  is  determined  by  the  relative  magnitudes  of  the 
driving  and  damping  mechanisms  within  the  combustor  that  add  and  remove  energy  from  the 
oscillations,  respectively.  If  the  acoustic  energy  gains  outweigh  the  losses,  small  disturbances 
amplify  to  levels  that  can  adversely  affect  guidance,  thrust  history,  and  the  structural  integrity  of 
the  motor'*^.  Because  of  the  high  energy  densities  and  low  losses  observed  in  most  combustors 
designed  for  use  in  propulsion  systems,  the  possibility  of  encountering  instability  is  high. 
Combustion  instability  is  a  common  problem  in  solid  propellant  rocket  motors^  as  well  as  in  other 
types  of  combustion  systems.  Since  combustion  instability  can  lead  to  system  and/or  mission 
failure,  the  development  of  reliable  methods  for  predicting  the  stability  of  solid  propellant  rocket 
motors  has  been  a  long  standing  goal  of  the  propulsion  community. 

In  order  to  eliminate  or  reduce  the  severity  of  acoustic  oscillations  resulting  from 
combustion  instabilities  in  solid  rocket  motors,  the  mechanisms  by  which  energy  is  added  to  and 
removed  from  these  oscillations  must  be  identified  and  understood.  The  interaction  of  the  solid 
propellant  combustion  process  with  the  acoustic  oscillations  in  the  combustor  is  generally 
accepted  as  the  source  of  the  energy  necessary  to  initiate  and  maintain  the  instability.  On  the 
other  hand,  processes  such  as  energy  transmitted  (out  of  the  motor)  through  the  choked  nozzle, 
viscous  dissipation,  and  flow  turning  tend  to  remove  energy  from  the  oscillations  and,  thus, 
stabilize  the  motor.  While  viscous  damping  is  generally  neglected  in  stability  analysis,  and  much  is 
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known  about  the  processes  of  nozzle  damping  and  combustion  driving,  relatively  little  is  known 
about  the  flow  turning  loss. 

The  velocity  of  combustion  products  leaving  the  burning  solid  propellant  is  generally  radial, 
as  shown  in  Fig.  1.  As  these  gases  move  away  from  the  propellant  and  are  entrained  in  the  core 
flow,  th^  are  turned  from  the  radial  direction  to  axial.  During  this  turning  process,  the 
combustion  products  acquire  axial  acoustic  energy  from  the  core  flow  oscillations.  This 
phenomenon  was  first  discussed  by  Culick^*"  who  argued  that  since  the  combustion  products 
entering  the  core  flow  acquire  acoustic  energy  from  the  existing  axial  acoustic  core  flow 
oscillations,  the  axial  acoustic  field  experiences  a  loss  of  acoustic  energy,  which  tends  to  stabilize 
the  system. 

The  flow  turning  phenomenon  has  been  studied  experimentally  by  Hersch"’’^  and  Chen*^ 
and  analytically  by  Baum'^’*^  and  Hegde  et  al'^‘*.  In  these  investigations,  the  behavior  of 
acoustic  waves  in  a  channel  in  which  the  solid  propellant  combustion  process  was  simulated  by 
mass  addition  through  a  side  wall  was  investigated.  Hersch"’*^  used  pressure  measurements 
upstream  and  downstream  of  a  region  of  wall  injection  to  determine  the  influence  of  flow  turning 
upon  induced  pressure  oscillations.  This  study  concluded  that  mass  addition  through  the  side  wall 
resulted  in  damping  of  the  axial  acoustic  field.  The  magnitude  of  the  flow  turning  loss  was 
determined  to  depend  upon  the  frequency  of  oscillation,  and  was  not  linearly  dependent  upon  the 
injection  rate.  These  conclusions  are  inconsistent  with  the  behavior  of  the  flow  turning  loss  as 
predicted  by  Culick^'°.  The  method  used  in  this  study  determines  the  overall  energy  loss  due  to 
the  injection,  but  cannot  isolate  the  flow  turning  loss  term  as  identified  by  Culick  from  other 
effects  of  the  same  order  of  magnitude  as  the  flow  turning  loss  that  were  present  in  the 
experiment.  Chen'^  measured  the  acoustic  intoisity  flux  through  the  boundaries  of  a  control 
volume  in  the  flow  turning  r^on  to  determine  the  damping  due  to  flow  turning.  This 
experimental  method  also  cannot  isolate  the  flow  turning  loss  from  other  effects  of  the  order  of 
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similar  order  that  are  coupled  with  it.  Baum's'^’*^  computational  studies  show  that  the  flow  field 
near  the  wall  is  extremely  complex,  involving  interactions  between  a  multitude  of  vortical 
structures  and  energy  exchange  between  the  mean  flow  and  the  oscillations.  Hegde'^'^  used  a 
linear  analysis  and  hot  wire  measurements  to  determine  the  acoustic  characteristics  and  losses  of  a 
duct  containing  a  standing  axial  acoustic  wave  and  mass  addition  through  the  side  wall.  This 
investigation  tested  the  validity  of  a  linear  stability  model  similar  to  those  currently  used  to  predict 
solid  propellant  rocket  stability  and  performance,  and  in  contrast  to  Baum's  results,  Hegde's 
studies  suggest  that  a  linear  model,  similar  to  those  currently  used  to  predict  the  stability  of  solid 
propellant  rocket  motors,  may  be  adequately  describe  the  important  features  of  the  flow  turning 
loss. 


The  following  questions  were  addressed  in  this  study; 

1 .  Is  the  flow  turning  loss  properly  included  in  currently  used  solid  propellant  rocket  motor 
stability  prediction  programs? 

2.  What  is  the  best  approach  for  measuring  the  flow  turning  loss? 

3.  How  does  the  flow  turning  loss  depend  upon  the  wall  injection  velocity,  the  axial  core  flow 
velocity,  the  fi-equency,  and  the  location  of  the  flow  turning  region  relative  to  the  standing 
pressure  wave? 

In  the  following  section,  the  accomplishments  of  the  research  performed  under  the  current 
program  are  discussed.  First,  the  theoretical  investigation  of  the  flow  turning  loss,  in  which  the 
basis  for  the  experimental  approach  was  developed,  is  described.  Second,  a  simple  computational 
model  used  to  study  the  behavior  of  the  flow  turning  loss  and  to  assist  in  the  interpretation  of 
experimental  results  is  discussed.  Finally,  the  re»jlts  of  the  experimental  study  are  presented. 
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RESEARCH  ACCOMPLISHMENTS 


1)  Theoretical  Investigation 

A  theoretical  investigation  of  the  flow  turning  loss  was  performed  in  order  to  determine 
whether  the  flow  turning  loss  is  accurately  included  in  state  of  the  art  stability  prediction  methods 
and  to  develop  a  practical  method  for  measuring  it.  First,  in  order  to  verify  the  existence  and 
mathematical  form  of  the  flow  turning  loss  as  defined  by  Culick^,  a  one-dimensional  acoustic 
stability  equation  was  derived  using  an  energy  balance  approach  similar  to  that  used  by  Cantrell 
and  Hart*^  and  later  by  Flandro^^  with  the  inclusion  of  flow  rotation  effects.  The  development  of 
the  stability  equation  begins  with  consideration  of  the  energy  equation  accurate  to  second  order  in 
acoustic  perturbations  and  first  order  in  Mach  number.  A  useful  expression  is  then  found 
containing  products  of  only  first  order  quantities  by  using  the  mass  and  momentum  equations. 
The  resulting  stability  equation  is  similar  to  that  originally  derived  by  Culick  using  a  different 
theoretical  approach.  The  stability  equation  for  a  two  dimensional  rectangular  volume  of  length  L 
and  height  /f  can  be  written  as 


Po<*o 


+ 


2T\  \  r^fcii 


H 

0 


where  a  is  the  acoustic  growth  rate  constant,  u  and  v  represent  the  velocity  components  in  the 
axial  and  transverse  directions  respectively,  and  m,,  is  the  mass  flux  at  the  wall.  A  numerical 
subscript  represents  the  order  of  perturbation  of  the  subscripted  term  (e.g.,  Pi  is  the  first  order 
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pressure  perturbation).  Angle  brackets,  <  >,  represent  an  average  over  the  subscript  variable 
(e.g.,  (  represents  a  time  average).  is  a  term  with  units  of  energy  given  by 


flo  2 


I  dx 


(2) 


Flandro^°  states  that  the  inclusion  of  the  second  term  in  Eq.  2  is  somewhat  controversial,  and 
demonstrates  that  this  term  is  correctly  placed  and  does  not  lead  to  results  that  are  inconsistent 
with  solutions  of  the  linear  problem  carried  out  by  other  means. 

The  organization  of  the  terms  of  Eq.  I  was  performed  in  a  way  that  facilitates  a  comparison 
with  Culick's  result.  Cuiick*^,  states  "the  combination  [of  terms]  can  be  written  in  many  equivalent 
ways;  it  is  by  no  means  obvious  at  this  point  which  is  most  suitable".  Culick  argues  that  his 
chosen  manner  of  combining  the  terms  aids  interpretation  and  incorporation  into  the  three- 
dimensional  problem.  The  term  identified  by  Culick  as  "mean  flow/acoustic  surface  interactions", 
later  to  become  known  as  flow  turning  term,  in  Eq.  1  is 


FLOW  TURNING  LOSS  = 


(3) 


Inspection  of  Eq.  3  shows  that  it  will  always  be  negative  in  a  rocket  motor,  where  mass 
enters  through  the  burning  surfaces.  For  such  a  configuration,  this  term  tends  to  reduce  a, 
thereby  tending  to  stabilize  the  system.  On  the  other  hand,  if  mass  exits  the  system  through  a  side 
wall  (such  as  through  a  T-bumer  nozzle),  the  term  is  positive,  which  increases  a  and  destabilizes 
the  system. 

This  analysis  confirmed  that  the  flow  turning  loss  term  is  appropriately  included  in  the  one 
dimensional  acoustic  stability  equation.  The  analysis  also  confirmed  Van  Moorhem's^'  conclusion 
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that  the  flow  turning  loss  is  not  absent  from  the  three  dimensional  analysis  as  suggested  by 
Culick^,  but  that  the  "classical  flow  turning"  term  can  only  be  isolated  after  cross  sectional 
averaging.  Tracing  the  origin  of  the  flow  turning  loss  term  backwards  through  the  derivation 
shows  that  it  originates  in  a  rotation  term  of  the  axial  momentum  equation.  The  flow  turning  loss 
is  implicitly  included  in  any  multi-dimensional  analysis  that  does  not  impose  irrotationality. 

The  one  dimensional  stability  equation  shown  in  Eq.  1  was  derived  by  expanding  the 
governing  equations  in  terms  of  first  and  second  order  quantities.  While  such  a  procedure  is 
theoretically  sound,  it  is  difficult  to  correlate  first  and  second  order  quantities  with  measured 
experimental  data.  Therefore,  a  new  theoretical  approach,  detailed  in  Ref  22,  has  been  developed 
which  is  more  useful  for  the  experimental  determination  of  the  flow  turning  loss.  This  approach 
starts  with  the  equations  of  mass  and  axial  momentum  conservation  in  which  each  of  the 
dependent  variables  is  expressed  in  terms  of  a  time  averaged  and  a  fluctuating  component.  For 
example,  the  pressure  is  expressed  as  p  =  p+p',  where  p  is  the  time  average  of  the  measured 
pressure  p. 

The  expression  obtained  by  this  method  for  the  complex  growth  rate  constant  of  the 
oscillations  is 


+  P 


(4) 


Comparison  of  Eqs.  1  and  4  reveals  that  the  form  of  the  equations  is  identical.  However,  the 
quantities  involved  in  these  equations  are  defined  differently.  Equation  3  lends  itself  more 
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rigorously  to  the  verification  of  experimental  results  due  to  the  way  in  which  the  terms  have  been 
defined. 


In  the  above  developments,  the  flow  was  assumed  isentropic  to  facilitate  the  use  of  the 
energy  balance  approach  based  upon  the  Cantrell  and  Hart  analysis.  The  result  of  this  assumption 
was  that  Eq.  4  is  limited  to  flows  with  uniform  mean  temperature  and  density.  The  above  results 
cannot  be  used  to  analyze  combustion  systems  in  which  mean  temperature  gradients  and  heat 
losses  may  be  significant.  In  order  to  determine  the  effects  of  temperature  gradients  upon  the 
acoustic  stability  and  the  flow  turning  loss  term,  an  investigation  of  the  non-isentropic  equations 
was  performed. 

In  this  analysis,  a  stability  equation  that  accounts  for  the  mean  temperature  and  density 
gradients  that  appear  in  the  hot  flow  experimental  studies  was  developed  from  the  conservation 
equations  for  mass,  axial  momentum,  and  energy: 


~  +  — (pi/)+— (pv)  =  0 


(5) 


du  du 

p_^p„_^pv 


du  dp 
— +— 
^  dx 


=  0 


(6) 


( ds  ds  cis^ 

(7) 

where  s  is  the  specific  entropy,  R  is  the  gas  constant  and  Q  is  the  rate  of  heat  addition  per 
volume.  The  energy  equation  was  not  used  explicitly  in  the  previous  development.  It  was 
implicitly  included  in  the  isentropic  relations.  These  equations  were  used  together  with  the 
perfect  gas  equation  of  state. 
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p  =  pRT 


(8) 


w^ere  T  *<’  the  temperature.  In  the  analysis,  it  was  assumed  that  Q'=0  (:.e.,  the  unsteady  heat 
addition  or  loss  within  the  volume  was  neglected).  The  development,  which  is  similar  to  the 
analysis  used  in  the  previous  section,  is  detailed  in  Ref.  23.  The  resulting  acoustic  stability 
equation  is; 


u'+puu'^ 


]t\  ti!  up  (  ^  i/pY 

/ tj,  'o  \  po^  1 j  p  i,  Sr  J 


dx 


where 


(9) 


(10) 


The  most  important  differences  to  note  between  £q.  9  and  Eq.  4  is  that  p  is  not  constant  and 
a*  ^  —  in  Eq.  9.  The  term 
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is  identified  as  the  flow  turning  term  in  this  version  of  the  stability  equation,  and  compares  to  the 
flow  turning  term  of  the  constant  density  stability  equation,  shown  in  Eq  .  3. 

The  term 

which  appears  in  Eq.  9,  does  not  appear  at  all  in  Eq.  4  because  [p]"  =  0  in  the  previous  analysis. 
Comparison  of  the  equations  shows  that  the  term 


in  the  stability  equation  that  accounts  for  temperature  gradients,  Eq.  9,  relates  to  the  term 


in  Eq.  4.  The  term  fi'om  Eq.  4  represents  fluxes  through  the  axial  boundaries  of  the  control 
volume,  while  the  term  from  Eq.  9  cannot  be  reduced  to  a  surface  flux  term,  and  represents 
effects  that  occur  inside  the  volume.  This  difference  is  due  to  the  inclusion  of  density  gradients  in 
Eq.  9,  and  suggests  that  the  interaction  of  axial  acoustics  with  axial  density  gradients  in  the 
presence  of  a  mean  flow  may  drive  or  damp  the  acoustics. 

2)  Numerical  Simulation 
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In  an  effort  to  determine  the  behavior  of  the  various  terms  involved  in  the  derived  acoustic 
stability  equation,  Eq.  3,  the  experimental  setup  used  in  the  flow  turning  investigation  was 
modeled  numerically.  Figure  2  shows  the  simulated  setup  in  which  the  length  L,  the  height  H,  the 
distance  from  the  head  end  to  the  burner  Zl,  and  the  distance  from  the  head  end  to  the  end  of  the 
burner  L2  are  characteristic  lengths  of  the  experimental  setup.  This  simulation  involved  the 
derivation  and  numerical  solution  of  the  mass  and  momentum  conservation  equations  that 
describe  the  experimental  setup  and  the  test  conditions.  For  this  analysis,  the  flow  was  assumed 
to  be  isentropic  with  perturbations  that  are  harmonic  in  time.  The  two  dimensional  forms  of  the 
equations  were  averaged  over  the  cross-section  of  the  duct  and  then  integrated  along  the  axis  of 
the  experimental  setup.  The  solutions  were  then  substituted  into  the  developed  acoustic  stability 
equation  to  determine  the  values  of  the  terms  that  affect  the  stability  of  the  system  and  to 
determine  whether  the  numerical  solution  satisfies  Eq.  4. 

It  should  be  noted  that  the  numerical  simulation  models  the  experimental  setup  and  not  an 
unstable  solid  propellant  rocket  motor.  Several  important  differences  exist  between  the  behaviors 
of  the  experimental  setup  and  an  actual  solid  rocket  motor  that  affect  the  nature  of  the  simulation. 
The  most  important  difference  is  that  the  experiment  is  an  externally  driven  system,  while  an 
unstable  rocket  motor  is  a  self  driving  system.  One  result  is  that  the  experimental  system  is 
always  neutrally  stable;  i.e.,  the  oscillations  neither  amplify  nor  decay  in  time.  This  means  that 
while  the  various  terms  of  the  acoustic  stability  equation  may  change  for  different  configurations, 
the  resulting  growth  rate  constant  a  must  equal  0  for  the  experiment.  On  the  other  hand,  the 
growth  rate  for  the  oscillations  in  a  rocket  motor  is  determined  by  the  sum  of  the  terms  in  Eq.  4. 
Another  difference  is  that  the  frequency  of  oscillation  in  the  experimental  setup  is  known;  i.e., 
chosen  by  the  experimenter,  whereas  the  frequency  of  oscillation  in  an  unstable  rocket  motor  is 
determined  by  relevant  conditions  (e  g.,  the  geometry  and  temperature)  and  is  not  easily 
predicted.  With  this  in  mind,  the  ability  of  the  numerical  simulation  to  satisfy  the  derived  stability 


equation  is  determined  by  summing  the  calculated  values  of  the  terms  on  the  right  hand  side  of 
Eq.  4.  The  sum  of  the  terms,  and  thus  the  growth  rate  a,  should  equal  0.  Calculations  were 
carried  out  for  a  number  of  different  experimental  configurations,  and  in  all  cases,  excellent 
agreement  between  the  results  of  the  numerical  simulation  and  the  developed  stability  equation 
was  observed. 

Results  of  a  simulation  for  a  simple  experimental  configuration  are  shown  in  Figs.  3  and  4. 
For  this  example,  the  admittances  of  the  upstream  injector,  the  top  and  bottom  walls,  and  the 
burner  surface  are  assumed  to  be  zero  in  order  to  simplify  the  calculations.  The  conditions  for  the 
simulations  have  been  set  as  follows;  1)  the  mean  pressure  is  atmospheric,  2)  the  amplitude  of  the 
pressure  oscillation  at  the  injector  face  is  200  Pa,  3)  the  burner  is  10  cm  long  and  is  located  70  cm 
downstream  of  the  injector,  and  4)  the  mean  Mach  number  of  the  flow  through  the  burner  is 
0.004,  which  simulates  the  Mach  number  of  the  flow  at  the  multi-difiusion  flame  burner  surface 
for  a  typical  cold  flow  experiment.  Figure  3  describes  the  axial  variations  of  the  amplitude  and 
phase  of  the  acoustic  pressure  oscillation  in  the  experimental  setup.  The  location  of  the  burner  is 
shown  in  the  figure.  It  is  seen  fi-om  the  figure  that  the  effect  of  mass  addition  upon  the  pressure  is 
small  for  the  investigated  Mach  number.  The  axial  variations  in  of  the  relevant  terms  of  the 
acoustic  stability  equation,  Eq.  4,  are  described  in  Fig.  4.  The  sum  of  these  terms,  which  is 
proportional  to  a,  is  also  shown  (term  6  in  the  figure),  and  is  equal  to  0  at  all  locations,  which  is 
in  agreement  with  the  theory  and  the  discussion  in  the  previous  paragraph.  It  should  be  noted  that 
for  this  burner  surface  Mach  number,  the  calculated  flow  turning  loss  is  small.  This  is  due  to  the 
low  Mach  number  and  the  short  burner  length  (the  flow  turning  loss  accumulates  over  length), 
and  does  not  suggest  that  the  flow  turning  loss  is  insignificant  in  a  solid  rocket  motor.  Figure  4 
also  shows  that  the  terms  of  Eq.  4  which  are  of  first  order  in  Mach  number  tend  to  be  of  similar 
magnitude.  As  a  result,  previous  experimental  studies'^  could  not  effectively  measure  the  flow 
turning  loss  because  the  effects  of  all  relevant  terms  were  not  accounted  for. 
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The  effects  of  several  parameters  upon  the  magnitude  of  the  flow  turning  loss  were  studied 
using  the  numerical  simulation  procedure.  One  such  parameter  was  the  location  of  the  simulated 
burner  relative  to  the  standing  acoustic  wave.  Figure  S  indicates  that  the  flow  turning  loss  is 
dependent  upon  the  relative  location  of  the  side  wall  mass  injection.  The  flow  turning  loss  is 
greatest  when  the  burner  is  located  at  an  axial  acoustic  velocity  maximum,  and  is  least  when  the 
burner  is  at  a  velocity  minimum.  Another  investigated  parameter  was  the  mean  core  flow  Mach 
number.  A  series  of  simulations  in  which  the  mean  core  flow  Mach  number  M^  was  varied  fi'om 
M^  =  0.0  to  M^  =  0.1  demonstrated  that  the  magnitude  of  the  flow  turning  loss  is  independent  of 
the  core  flow  velocity  over  the  investigated  range.  Figure  6  shows  that  while  the  magnitude  of 
some  terms  change  (specifically,  the  mean  flow  convection  terms),  the  flow  turning  term  remains 
unchanged  in  all  cases.  Results  of  the  third  parameter  studied,  the  burner  injection  Mach  number 
Mg,  are  displayed  in  Fig.  7.  This  investigation  showed  that  the  magnitude  of  the  flow  turning  loss 
varies  linearly  with  the  injection  Mach  number. 

The  conditions  used  for  the  iiutial  numerical  simulations  were  idealized  in  order  to  simplify 
the  interpretation  of  the  results.  The  actual  conditions  encountered  in  the  experimental 
investigation  are  rather  complex,  and  difficult  to  approximate  in  a  numerical  simulation. 
Therefore,  the  numerical  simulation  was  modified  to  use  experimentally  determined  (i.e., 
measured)  boundary  conditions.  This  version  of  the  code  is  used  to  aid  in  interpretation  of 
experimental  results. 

The  simulations  involving  experimentally  determined  boundary  conditions  require  measured 
values  of  the  mean  axial  velocity  component  at  the  initial  axial  boundary  of  the  control  volume  as 
shown  in  Fig.  8,  the  amplitude  and  phase  of  the  axial  acoustic  velocity  component  at  both  axial 
boundaries,  as  well  as  the  vertical  velocity  component  and  the  amplitude  and  phase  of  the  vertical 
acoustic  velocity  along  the  top  and  bottom  boundaries.  A  shooting  method  is  then  used  to 
determine  the  acoustic  velocity  and  pressure  solutions  which  minimize  the  errors  at  the  axial 
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boundaries  in  a  least  squares  sense.  These  values  of  velocity  and  pressure  are  then  used  along 
with  Eq.  3  to  predict  values  of  the  various  terins  in  the  acoustic  stability  equation.  The  results  of 
this  type  of  simulation  and  the  comparison  with  experimental  results  are  discussed  in  the  next 
section. 

3)  Experimental  Study 

The  experimental  setup  used  in  this  study  is  shown  schematically  in  Fig.  9.  Details  of  the 
experimental  setup  are  given  if  Ref  24.  The  setup  is  designed  to  simulate  flow  phenomena  near 
the  walls  of  a  solid  propellant  rocket  experiencing  an  axial  acoustic  instability.  In  the  first  phase 
of  the  study,  cold  flow  experiments  were  performed  to  determine  the  validity  of  the  "flow  turning" 
loss  and  its  effect  upon  the  acoustic  stability.  In  the  second  phase  of  the  experimental  study,  a 
multi-diffusion  flame  burner  (shown  schematically  in  Fig.  10)  was  used  to  better  simulate  actual 
motor  conditions  and  to  determine  the  influence  of  temperature  gradients  upon  the  flow  turning 
loss  and  upon  other  factors  relevant  to  the  stability  of  rocket  motors. 

a)  Cold  Flow  Studies 

In  the  cold  flow  investigation,  the  flow  turning  loss  and  were  determined  by  measuring 
the  two  components  of  the  mean  and  acoustic  velocities  and  the  acoustic  pressure  and  then 
substituting  the  measured  values  into  Eq.  4.  The  term  describes  the  contribution  of  the  flow 
turning  loss  to  the  growth  rate  constant  a  of  the  investigated  re^on,  and  is  defined  as: 

>^Ko]o 

A  number  of  experiments  were  performed  for  various  conditions  to  determine  the  effects  of  these 
conditions  upon  the  flow  turning  loss.  Several  hours  were  required  to  acquire  the  necessary 
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velocity  and  pressure  data  for  each  configuration.  Due  to  the  slow  data  acquisition  and  the 
sensitivity  of  the  flow  field  to  experimental  conditions,  only  a  small  number  of  data  per 
configuration  (generally  3)  were  collected.  To  compensate,  a  large  number  of  experimental  runs 
were  performed,  and  trends  were  determined  by  comparing  a  number  of  results. 

The  results  displayed  in  Fig.  1 1  shows  the  effect  of  the  burner  injection  velocity  upon 
in  the  investigated  region.  The  results  were  obtained  in  a  test  conducted  with  a  fi-equency  of  550 
Hz.,  a  mean  core  flow  velocity  of  0.4  m/s,  and  injection  at  a  pressure  node.  The  measured 
varies  linearly  with  the  mean  burner  injection  velocity. 

Figure  12  shows  the  dependence  of  in  the  investigated  region  upon  the  mean  core  flow 
velocity.  This  test  was  conducted  with  a  fi’equency  of  550  Hz.,  a  burner  injection  velocity  of 
0. 156  m/s,  and  the  burner  located  at  a  pressure  node.  The  measured  increase  with  the  mean 
core  flow  velocity.  The  trend  of  this  data  and  the  results  of  other  tests  (because  a  conclusion 
cannot  be  drawn  from  three  data  points)  shows  that  appears  to  approach  a  limit  value  as  the 
core  velocity  increases.  This  is  due  to  incomplete  turning  of  the  flow  inside  the  investigated 
region.  In  the  experimental  study,  the  investigated  control  volume  does  not  extend  from  the 
burner  surface  to  the  top  wall.  Inspection  of  the  mean  velocity  profile  shown  in  Fig.  13  reveals 
that  the  flow  at  the  top  of  the  investigated  control  volume  is  not  fully  turned  into  the  axial 
direction.  As  the  core  flow  Mach  number  increases,  flow  turning  is  completed  to  a  higher  degree 
within  the  investigated  control  volume.  Therefore,  it  is  expected  that  as  the  core  flow  Mach 
number  increases,  the  flow  turning  loss  will  also  increase,  which  is  consistent  with  the  trend  of  the 
measured  data. 

The  effect  of  varying  the  amplitude  of  the  excited  standing  wave  upon  is  displayed  in 
Fig.  14.  The  results  are  fi-om  an  experiment  in  which  the  mean  core  flow  velocity  was  1.45  m/s, 
the  mean  burner  injection  velocity  was  0.62  m/s,  the  frequency  was  550  Hz.,  and  the  injection 
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occurred  at  an  acoustic  pressure  antinode.  The  pressure  amplitudes  were  normalized  with  respect 
to  the  maximum  measured  amplitude.  The  f.gure  shows  that  is  independent  of  the  amplitude 
of  oscillation,  which  is  expected  in  a  linear  analysis. 

The  dependence  of  upon  the  location  of  the  flow  turning  region  with  respect  to  the 
standing  acoustic  wave  is  shown  in  Fig.  15.  These  data  were  obtained  in  experiments  with  a  mean 
core  velocity  of  2.86  m/s,  a  mean  burner  injection  velocity  of  1.34  m/s,  and  a  fi’equency  of  550 
Hz.  The  measured  was  maximum  when  the  flow  turning  occurred  at  an  acoustic  pressure 
node,  and  minimum  when  the  flow  turning  region  was  located  at  a  pressure  antinode. 

b)  Comparison  of  Experimental  Data  to  Calculations 

As  discussed  above,  the  numerical  scheme  was  modified  to  use  measured  boundary 
conditions  to  aid  in  the  interpretation  of  the  measured  data.  Measured  and  calculated  values  of 
the  axial  acoustic  velocity  for  an  experimental  test  are  shown  in  Fig.  16.  The  agreement  between 
the  measured  and  calculated  velocity  values  is  excellent.  Figure  17  shows  comparisons  between 
measured  and  calculated  acoustic  pressures  and  measured  and  calculated  phase  angle  between  the 
pressure  and  velocity  oscillations.  A  good  qualitative  agreement  between  the  calculated  and 
measured  results  is  observed,  but  the  error  in  the  predicted  phase  values  is  large.  The  calculated 
and  measured  values  of  the  flow  turning  loss  for  this  test  configuration  are  shown  in  Fig.  18.  An 
excellent  agreement  is  observed  between  the  measured  and  calculated  values.  Note  that,  for  the 
experimental  configuration  whose  results  are  shown  in  this  section,  the  value  of  the  flow  turning 
loss  term  was  negative,  indicating  a  gain  of  acoustic  energy.  This  is  the  case  because  of  the 
nature  of  the  flow  in  the  investigated  control  volume.  The  flow  injected  through  the  burner  is 
turned  into  the  axial  direction,  but  at  the  same  time,  the  core  flow  entering  the  control  volume 
fi’om  upstream  is  being  turned  upward.  For  some  experimental  configurations,  the  flow  turning 
loss  due  to  the  turning  of  the  injected  flow  is  smaller  than  the  'flow  turning  gain'  due  to  the  turning 
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of  the  core  flow.  This  flow  turning  gain  effect  is  observed  only  because  of  the  choice  of  control 
volumes;  in  a  volume  extending  fi-om  wall  to  wall  (e.g.,  in  a  solid  rocket  motor)  flow  turning  will 
always  result  in  a  loss  of  acoustic  energy. 

c)  Hot  Flow  Studies 

The  results  of  an  experimental  study  of  the  effects  of  mean  temperature  gradients  upon  the 
flow  turning  loss  and  other  factors  relevant  to  the  stability  of  rocket  motors  are  discussed  in  this 
section.  Combustion  is  present  in  these  hot  flow  studies,  but  the  combustion  region  itself  is  kept 
external  to  the  investigated  region,  as  described  in  Fig.  19.  The  driving  and  damping  effects  of 
the  combustion  process  have  been  studied  previously  in  this  research  program  by  Sankar^^  and 
Chen’^,  and  the  inclusion  of  these  mechanisms  would  serve  to  complicate  the  flow  turning  analysis 
by  coupling  additional  terms  in  the  stability  analysis.  The  data  acquired  in  the  hot  flow 
experiments  was  analyzed  using  the  form  of  the  stability  equation  shown  in  Eq.  9,  which  accounts 
for  spatial  variations  in  the  mean  density. 

A  large  temperature  difference  exists  between  the  core  flow  (which  consists  of  room 
temperature  air)  and  the  flow  of  hot  combustion  products  leaving  the  flat  flame  that  is  stabilized 
above  the  burner  surface.  This  generates  a  line  of  maximum  temperature  gradient  begins  at  the 
upstream  edge  of  the  burner  and  follows  the  turning  of  the  injected  flow,  as  illustrated  in  Fig.  20. 
The  flow  turning  loss  and  in  the  hot  flow  experiments  were  determined  by  measuring  the  two 
components  of  the  mean  and  acoustic  velocities,  the  mean  density,  and  the  acoustic  pressure  in 
the  investigated  region  and  then  substituting  the  measured  values  into  Eq.  9.  Experiments  were 
performed  for  various  conditions  to  determine  the  effects  upon  the  flow  turning  loss. 

The  results  shown  in  Fig.  21  were  obtained  in  a  test  conducted  with  a  fi-equency  of  580  Hz., 
a  mean  core  mass  flow  rate  of  3.63  g/s,  and  the  burner  located  at  a  pressure  node,  and  cold  flow 
results  fi’om  an  experiment  in  which  the  core  mass  flow  rate  was  the  same,  the  burner  was  located 
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at  a  pressure  node,  and  the  frequency  was  reduced  to  300  Hz.  The  frequency  was  reduced  in  the 
cold  flow  tests  to  preserve  the  wavelength  of  the  pressure  oscillation  under  the  cold  flow 
conditions.  The  linear  dependence  of  upon  the  burner  mass  flow  rate  agrees  with  the  theory. 
Comparison  of  the  hot  and  cold  results  shows  that  the  overall  flow  turning  loss  is  greater  in  the 
cold  flow  cases.  This  is  somewhat  misleading,  however,  because  the  geometry  of  the  hot  flow 
Add  is  not  the  same  as  that  of  the  cold  flow  field.  While  the  mass  flow  rates  are  preserved 
between  the  two  corifigurations,  the  velocity  out  of  the  burner  surface  is  much  higher  in  the  hot 
flow  case.  The  result  of  the  flow  differences  is  that  the  stream  of  hot  combustion  products 
leaving  the  burner  turns  more  gradually  as  it  moves  away  from  the  burner  than  the  flow  of 
injected  cold  flow  air.  A  qualitative  comparison  can  be  made,  however,  between  the  behavior  of 

under  hot  and  cold  conditions  to  determine  the  effect  of  the  presence  of  temperature 
gradients.  Qualitatively,  the  results  show  that  the  effect  of  the  wall  injection  velocity  upon  flow 
turning  loss,  as  determined  by  is  virtually  unaffected  by  the  presence  of  mean  temperature 
and  density  gradients  in  the  flow  turning  region.  The  theoretical  analysis  shows  that  some  effect 
of  the  density  gradient  should  be  present,  but  the  measurements  suggest  that  the  influence  of  the 
density  gradient  is  small  compared  to  the  effect  of  the  injected  mass  flow  rate  in  the  studied 
configuration.  It  is  likely  that  the  effect  of  the  density  gradients  upon  the  flow  turning  loss  in 
actual  solid  rocket  motors  >^11  also  be  small,  because  mean  density  changes  encountered  in  actual 
solid  rocket  motors  are  generally  small  fi-actions  of  the  total  density. 

Figure  22  shows  the  dependence  of  upon  the  mean  core  mass  flow  rate  measured  in 
tests  conducted  with  a  fi-equency  of  580  Hz.  and  an  injected  mass  flow  rate  of  2.26  g/s  through 
the  burner,  which  was  located  at  a  pressure  node.  The  measured  increase  with  the  mean 
core  flow  velocity  and  appears  to  approach  a  limit.  This  behavior  results  from  incomplete  flow 
turning  in  the  investigated  region,  as  was  explained  in  the  discussion  of  the  cold  flow  results.  The 
cold  flow  data  is  fi’om  an  experiment  in  which  the  mass  flow  rate  through  the  burner  was  the 
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same,  the  burner  was  located  at  a  pressure  node,  and  the  frequency  was  reduced  to  300  Hz  (to 
preserve  the  wavelength).  Qualitative  comparison  again  shows  little  difference  between  the 
behavior  of  in  the  hot  and  cold  flow  experiments. 

The  dependence  of  a_  upon  the  location  of  the  flow  turning  region  with  respect  to  the 
standing  acoustic  wave  is  shown  in  Fig.  23.  These  data  were  obtained  in  hot  and  cold 
experiments  with  a  mean  core  mass  flow  of  3.63  g/s,  and  a  mean  mass  flow  through  the  burner  of 
2.26  g/s.  The  measured  was  maximum  when  the  burner  was  located  at  an  acoustic  pressure 
node,  and  minimum  when  the  burner  was  located  at  a  pressure  antinode.  No  influence  of  the 
density  gradient  upon  the  behavior  of  is  suggested  by  the  theory  in  this  case,  because  the 
position  of  the  burner  relative  to  the  standing  acoustic  wave  has  no  influence  on  the  mean 
temperature  distribution  in  the  combustor. 

The  terms  of  the  acoustic  stability  equation  that  are  proportional  to  the  mean  Mach  number 
(e.g.,  flow  as  turning  loss  term)  measured  in  this  study  were  found  to  be  much  smaller  than  the 
Mach  number  independent  terms  (i.e.,  the  classical  acoustic  intensity  terms)  in  the  experimental 
setup.  This  is  due  to  the  large  acoustic  admittance  of  the  burner,  which  absorbs  a  large  fraction  of 
the  acoustic  energy  of  the  investigated  control  volume.  This,  however,  is  not  typical  of  the 
conditions  encountered  in  actual  solid  propellant  rocket  motors,  in  which  the  response  of  the 
propellant  is  much  smaller  than  that  of  the  burner  used  in  this  study.  The  Naval  Weapons  Center 
has  been  involved  in  a  program  to  develop  an  increased  understanding  of  combustion  instability 
through  which  a  database  of  motor  and  stability  data  is  being  coUected^^^°.  Blomshield  et  al^° 
test  fired  23  motors  in  a  number  of  configurations,  measured  the  decay  rate  of  induced 
oscillations,  and  compared  the  data  to  stability  predictions  made  using  the  Standard  Stability 
Prediction  Program  (SSP-ID)^V  The  stability  analysis  performed  in  the  SSP  module  is  based 
upon  Culick's  stability  formulation.  The  results  show  that  for  the  test  motors,  the  flow  turning 
loss  is  very  important  to  the  overall  stability  of  the  motor.  Other  calculations  of  the  stability  of 
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full  scale  motors^'  show  that,  unlike  the  laboratory  combustor  used  in  this  study,  the  flow  turning 
loss  is  often  a  determining  factor  in  the  stability  of  solid  propellant  rocket  motors  due  to  the 
relatively  small  acoustic  response  of  solid  propellant  compared  to  the  multi-diffusion  flame  burner 
used  in  this  study. 

d)  Effect  of  Temperature  Gradients  Upon  The  Acoustic  Intensity 

In  a  study  performed  by  Chen*^,  a  significant  amount  of  acoustic  energy  was  absorbed  in  a 
region  above  a  side  wall  stabilized,  fiat  difilision  flame,  but  no  significant  energy  loss  was 
observed  under  cold  flow  conditions.  The  greater  damping  in  the  presence  of  combustion 
suggested  that  some  process,  possibly  due  to  the  presence  of  mean  temperature  gradients, 
contributes  to  the  dissipation  of  acoustic  energy  in  the  system.  This  result  does  not  agree  with 
analytical  studies,  which  show  that  for  low  mean  Mach  numbers,  the  presence  of  a  mean 
temperature  gradient  neither  drives  or  damps  acoustic  energy. 

An  experimental  study  of  Chen's  result  was  performed.  The  oscillatory  velocities  and 
pressures  were  measured  on  the  border  of  a  control  volume  located  above  the  burner  surface  as 
shown  in  Fig  24.  These  measurements  were  used  to  calculate  the  classical  acoustic  intensity  flux, 
entering  and  leaving  the  boundaries  of  a  two  dimensional  control  volume  above  the 

multi-diffusion  flame  burner.  The  control  volume  investigated  for  this  study  was  a  6.8  x  2.24  cm^ 
rectangle  located  at  the  center  plane  of  the  duct.  Measurements  of  the  oscillatory  velocities  and 
pressures  in  this  study  were  carried  out  with  the  burner  at  various  locations  relative  to  the 
standing  acoustic  wave,  with  and  without  combustion.  The  acoustic  oscillation  was  driven  at  700 
Hz. 

The  measured  acoustic  intensity  fluxes  were  used  to  determine  the  influence  of  the  location 
of  injection  relative  to  the  standing  acoustic  wave  upon  the  energy  absorption  in  the  control 
volume.  Figure  25  compares  the  ratio  of  the  total  acoustic  intensity  flux  out  of  the  control 
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volume  to  the  flux  through  the  downstream  boundary  (the  direction  of  the  acoustic  driving)  for 
hot  and  cold  flow  cases.  Figure  25  does  not  show  substantially  more  energy  absorption  in  hot 
flow  cases  than  in  cold  flow,  and  does  not  show  a  definite  dependence  of  the  energy  absorption 
upon  the  location  of  the  burner  in  the  acoustic  field.  An  average  error  of  one  degree  in  the  phase 
angle  between  the  pressure  and  velocity  can  result  in  as  much  as  a  S  percent  error  in  the 
calculated  energy  loss.  Within  this  range,  there  is  no  discernible  difference  between  the  energy 
absorption  with  and  without  combustion,  which  agrees  with  the  theoretical  results. 
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SUMMARY 


An  investigation  of  the  processes  that  control  the  flow  turning  loss  and  its  role  in  the 
damping  of  combustion  instabilities  in  solid  propellant  rocket  motors  has  been  performed. 
Theoretical,  experimental,  and  numerical  approaches  have  been  used  to  study  the  behavior  of  the 
flow  turning  loss  in  a  laboratory  combustor  designed  to  simulate  important  features  of  a  solid 
rocket  motor  experiencing  axial  instability. 

Theoretical  an^ysis  performed  in  this  study  has  confirmed  that  the  flow  turning  loss  term 
must  be  included  in  the  one  dimensional  acoustic  stability  equation.  The  analysis  also  confirmed 
that  the  flow  turning  loss  is  not  absent  fi'om  multi-dimensional  formulations,  but  that  the  term 
traditionally  identified  as  the  flow  turning  loss  occurs  as  a  consequence  of  averaging  the  equations 
over  the  cross  section  of  the  chamber.  The  spatial  avera^g  process  has  a  tendency  to  make  the 
physical  interpretation  of  terms  somewhat  difficult.  Tracing  the  flow  turning  loss  term  backwards 
through  the  analysis  has  shown  that  it  originates  in  the  rotational  terms  of  the  conservation 
equations,  and  indeed,  the  term  does  not  appear  in  analyses  that  do  not  account  for  rotation  of  the 
flow. 

An  acoustic  stability  equation  that  accounts  for  the  presence  of  mean  temperature  and 
density  gradients  in  the  chamber  has  been  developed.  This  equation  is  similar  in  form  to  the 
stability  formulae  developed  under  the  previous  analyses.  An  additional  term  appears  in  the 
stability  equation  that  accounts  for  cross-duct  temperature  gradients. 

An  experimental  study  was  conducted  in  a  setup  designed  to  simulate  flow  phenomena  near 
the  walls  of  a  solid  propellant  rocket  experiencing  an  axial  acoustic  instability.  A  cold  flow 
investigation  was  performed,  in  which  the  flow  turning  loss  and  were  determined  by 
measuring  the  two  components  of  the  mean  and  acoustic  velocities  and  the  acoustic  pressure  and 
then  substituting  the  measured  values  into  the  acoustic  stability  equation.  Experiments  were 
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performed  to  determine  the  effects  of  various  parameters  upon  the  flow  turning  loss.  The  flow 
turning  loss  was  shown  to  be  linearly  dependent  upon  the  wall  injection  velocity,  which  relates  to 
the  burning  rate  of  the  propellant  on  the  walls  of  a  rocket  motor.  This  result  agrees  with 
theoretical  predictions.  Measurements  showed  that  the  flow  turning  loss  in  the  experimental 
setup  increased  towards  a  limiting  value  as  the  core  flow  velocity  was  increased.  This  does  not 
agree  with  theoretical  predictions,  which  indicate  that  the  flow  turning  loss  is  independent  of  the 
core  flow  velocity.  This  discrepancy  is  due  to  incomplete  turning  of  the  injected  flow  in  the 
investigated  region,  and  that  under  the  conditions  encountered  in  an  actual  rocket  motor,  the  flow 
turning  loss  will  be  independent  of  the  core  flow  velocity.  The  flow  turning  loss  has  been  shown 
to  be  independent  of  the  amplitude  of  oscillation  over  the  investigated  range.  Because  terms  of 
third  order  and  higher  in  the  energy  equation  have  been  neglected,  this  behavior  cannot  be 
extrapolated  to  large  amplitude  oscillations.  This  is  of  no  real  concern,  however,  because  linear 
analysis  is  generally  performed  when  the  oscillations  are  small.  Also,  the  study  has  shown  that  the 
flow  turning  loss  is  maximum  and  minimum  when  the  side  wall  injection  occurs  at  acoustic 
pressure  nodes  and  antinodes,  respectively. 

The  terms  proportional  to  the  mean  Mach  number  (such  as  the  flow  turning  loss  term) 
measured  in  this  analysis  were  found  to  be  much  smaller  than  the  Mach  number  independent  terms 
(the  classical  acoustic  intensity  terms)  in  the  experimental  setup.  This  is  because  the  burner  is 
acoustically  'soft'  and  absorbs  a  great  deal  of  acoustic  energy.  The  result  is  that  the  zeroeth  order 
(Mach  number  independent)  losses  in  the  experimental  setup  are  high.  This  is  not  necessarily  the 
case  in  actual  solid  propellant  rocket  motors,  where  the  propellant  response  is  generally  small  and 
the  Mach  numbers  are  somewhat  higher.  In  actual  motors,  the  flow  turning  loss  and  other  terms 
of  similar  order  are  often  of  the  same  magnitude  as  the  propellant  response  terms. 

In  the  hot  flow  study,  the  effects  of  mean  temperature  gradients  upon  the  flow  turning  loss 
and  other  factors  relevant  to  the  stability  of  rocket  motors  were  investigated.  The  data  acquired 
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in  the  hot  flow  experiments  were  analyzed  using  the  stability  equation  formulated  to  account  for 
spatial  variations  in  the  mean  temperature.  Measurements  show  that  the  behavior  of  the  flow 
turning  loss  in  the  experimental  setup  is  not  significantly  affected  by  the  presence  of  a  mean 
temperature  gradient. 

In  contrast  to  the  results  of  the  previous  study,  measurements  did  not  show  substantial 
differences  in  the  amount  of  energy  absorbed  in  hot  flow  cases  to  that  absorbed  in  cold  flow,  and 
does  not  show  a  definite  dependence  of  the  energy  absorption  upon  the  location  of  the  burner  in 
the  acoustic  field.  The  current  results  show  that  there  is  no  difference  outside  the  range  of 
experimental  error  between  the  amount  of  energy  absorbed  with  and  without  combustion.  Theory 
suggests  that,  when  terms  proportional  to  the  mean  velocity  are  neglected,  the  stability  of  a  duct 
is  not  influenced  by  the  presence  of  mean  temperature  gradients,  which  agrees  with  the  results  of 
this  study. 
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ACOUSTIC  DRIVER 


Figure!.  Schematic  oftheconfijuntion  used  in  the  numerical  emulation  of  the 
experimental  setup. 
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Figure  3.  Numerically  predicted  amplitude  and  phase  of  the  pressure  oscillation  with  side 
ngection  near  an  acoustic  pressure  minimum. 
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Figure  4.  Numerical  prediction  of  the  relevant  tenns  of  the  acoustic  stability  equation. 
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Kgure  5.  Predicted  flow  turning  loss  for  injection  tt  en  acoustic  vdochy  maximum  and 
minimum,  with  the  acoustic  vdodty  shown  for  reference. 
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Figure  6.  Flow  turning  loss  and  remaining  terms  of  the  acoustic  stability  equation  for  3 
with  varying  core  flow  Mach  numbers. 
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Figure  9.  A  schematic  of  the  experimental  setup  utilized  in  the  investigation  of  flow 
turning  losses. 
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Figure  11.  Results  of  an  experiment  ihowing  the  effect  of  the  tide  wall  injection  velocity 
upona^. 


Figure  12.  Resuhs  of  an  eiqieriinem  showing  the  effect  of  the  mean  cm  flow  vdodty 
uponttj^. 
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Figure  13.  Vector  plot  of  t  typical  measured  mean  vdocHy  field  in  the  flow  turning 
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Figure  IS.  RMuht  of  an  eiqieriment  diowing  the  effect  ofthe  loctdon  of  the  flow  turniog 
region  relative  to  the  itanding  acoustic  wave  upon 


Figure  16.  Conipariwnoftlie  measured  uidcalcutatedfmgnitudes  of  the  anal 
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Figure  19.  A  schenuiticofthe  experimental  letup  and  the  location  of  tlie 
inveitigated  r^on. 
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Kgure24.  MeutumeMs  used  for  determiniiig  the  icoustic  intensity  entering  and  leaving 
foe  investigated  region. 
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figure  25.  Comparison  of  the  ratio  of  the  total  acoustic  intensity  flux  out  of  the  control 
vohime  to  the  flux  through  the  downstream  boundary  for  hot  and  cold  flow  cases. 
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field  in  n  anstable  aolid  rodcet  motor.  Although  die 
flow  turning  loss  is  ineotporated  in  sute  of  die  art 
stabiliqr  codes  of  solid  ncket  motots.  the  aaistenos  of 
this  lost  hat  never  been  adequately  verified 
experimentally.  Ihit  p^ier  describes  a  theoreUcal 
investigation  that  hat  been  used  to  guide  the 
development  of  an  experimental  approach  for  the 
determination  of  the  flow  turning  lots.  First,  the 
existence  of  the  flow  taming  lou  term  in  die  one- 
dimensional  subility  aquation  was  verified  using  an 
independent  energy  balance  qiproach.  Next,  a 
theoretical  model  that  could  be  used  to  determine  the  flow 
turning  loot  from  measured  dau  was  developed.  An 
equation  describing  the  acoustic  subility  of  the 
experimental  setup  in  terms  of  measurable  time  averaged 
and  fluctuating  quantities  was  derived  from  the 
conservation  equations.  Finally,  the  governing 
equations  for  die  experimental  setup  were  numerically 
integrated  and  the  eolutiont  substituted  into  the 
developed  ecouttic  subility  aquation  to  determine  die 
behavior  of  the  various  terms  that  •ppur  in  the  equation 
and  whether  the  numerical  solution  satisfies  die 
developed  expression  for  the  growdi  rau  of  the 
oscillations  in  the  experiment.  Calculations  were  carried 
out  for  a  number  of  different  experimental  configurations 
and  operating  conditions  which  showed  excellent 
agreement  between  the  resulu  of  the  numerical 
simulation  and  the  developed  stability  equation. 

Introduction 

In  Older  to  dimbiau  or  reduce  the  the  occunence  of 
combustion  instabilities  in  solid  propellant  rocket 
motors,  the  processes  diet  edd  or  remove  energy  from  die 
oscillatioiu  must  be  accurauly  accounted  for  in  a 
subility  analysis.  To  attain  this  goal,  an  understanding 
of  these  processes  must  be  developed.  It  is  generally 
accepted  that  the  energy  required  to  biitiaie  and  maintain 
the  insubility  is  supplied  by  the  combustion  process.  It 
is  elso  recognized  that  nozzle  damping,  viscous 
dissipation,  heat  transfer,  and  "flow  tuming*  are 
processes  dut  remove  energy  from  the  acoustic  field. 
This  study  is  concerned  with  die  understanding  and 
quantification  of  the  damping  provided  by  thv  "Oo* 
Aiming”  procen. 

~  Flow  turning  losses  in  an  nnsuble  rocket  motor  are 
caused  by  the  interaction  of  the  radially  directed  flow  ol 
combustion  products  that  leave  the  burning  propellant 
surface  with  the  axial  acoustic  field.  The  existence  of 


diis  flow  Biming  loss  wu  first  predicted  by  Culick^. 
Although  this  flow  turning  loss  is  curreatiy  iaootporaied 
in  stale  of  the  art  sUbiliqf  models  that  predict  the 
perfoimanoe  of  solid  propellam  rockeu,  itt  existence 
has  never  been  quantiutively  verified.  Previous 

affotts^^  under  diis  program  and  dsewhers  to  sMasure 
die  flow  taming  loes  have  failed  because  the  quantities 
that  needed  to  be  measured  were  of  die  same  order  as  the 
accuracy  of  the  measurement  tadmiques.  Consaquendy, 
this  diaotetieal  study  hat  been  undotaken  to  dmeimine 
whedier  die  .previously  developed  axpraesion  for  the 
flow  taming  lots  is  aecuraie  and  vdiether  dm  previous 
axperimental  studies  failed  to  account  for  all  the 
necessary  terms  in  their  evaluation  of  the  flow  turning 
lees,  ths^ore  giving  amneout  resulu. 

In  die  initial  phase  of  diis  program  an  attempt  was 
made  to  measure  the  flow  turning  lou  within  a  set  of 
eontrol  volumu  where  flow  turning  occuned.  This 
procedure  required  that  an  acoustic  energy  balance  be 
performed  on  die  investigated  control  volume.  In  this 
cue.  the  necessary  acoustic  energy  fluxes  were 
deiemiinad  from  measurementt  of  the  clauical  acoustic 
energy  flux  p‘*u',  vdiere  p'  and  n*  are  the  acoustic  pressure 
and  velociqr,  respectively,  acrou  all  the  boundariu  of 
die  bivutigued  control  volume.  Since  no  other  sources 
or  sinks  of  acoustic  energy  were  pruent  within  the 
faivutigatcd  control  volumu,  any  net  loss  of  acoustic 
anergy  wu  attributed  to  flow  tuming.  Unfortunately,  the 
lossu  measured  by  diis  procedure  were  of  die  same  order 
u  the  experiment^  enors  and.  consaquendy,  could  not 
be  Busted.  It  should  be  also  noted  that  a  ralai^  study  by 
Hcrsdi  and  Tso^  also  fafled  to  provide  reliable  rasuhs  due 
to  experimenul  biaccuraciu.  The  theoretical  study 
reported  herein  wu  stimulated  by  die  failure  of  all  studiu 
to  date  ta  accurately  measure  the  flow  tuming  lou. 

The  purpose  of  this  invutigation  is  to  determine 
the  causu  of  Qiea's^  foilure  to  measure  die  flow  turning 
lou  and  to  propow  an  alternate  experimental  approach. 
First,  in  erdu  to  verify  the  existence  of  the  flow  tuming 
lou.  a  one-dimensional  acoustic  suUlity  equation  wu 
developed  using  an  anergy  balance  qipioach  similar  to 
that  used  by  Cantrell  and  Han!^.  The  analysis  wu  based 
on  explicit  considerations  of  the  mau,  momentum,  and 
acoustic  energy  fluxu  at  the  control  volume  boundariu. 
It  involved  consideration  of  the  first  and  second  order 
perturbations  of  the  acoustic  quantities.  The 
considuation  of  second  order  quantitiu  is  neceuaiy  for 
proper  evaluation  of  terms  of  the  ordu  of  the  acoustic 
energy.  The  anafysis  started  with  an  expansion  of  the 
anergy  aquation  to  second  order  in  the  acoustic 
pertiulMtions.  Next,  second  order  expansions  of  the 
mus  and  momentum  conservation  aquations  were  used  to 
modify  die  anegy  aquation.  This  analysis  yielded  an 
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•xpressioii  for  the  growth  rate  a  of  die  omuble  wave 
that  if  independent  of  any  tcnns  involving  second  order 
acoustic  quantities.  The  resulting  stability  aouation  is 

similar  to  that  originally  derived  by  Culkk^  nsing  a 
different  dworedcal  approach. 

Subsequently,  in  order  to  develop  a  theoretical 
asodel  that  could  be  used  to  guide  die  development  of  an 
alternate  experimental  approach  for  determining  the  flow 
tumiiig  lou,  an  analysis  was  developed  that  allows  for 
the  expansion  of  all  dependent  variables  (e.g..  the 
pressure)  in  terms  of  time  averaged  and  fluctuating 
quantities  similar  to  those  measured  in  related 
experiments.  For  example,  the  preuure  is  described  by 
the  following  expression: 
p  -  p  +  p* 

where  P  is  the  time  average  of  die  measured  pressure,  P. 
These  dependent  variables  were  substituted  into  the  mus 
and  momentum  conservation  equatioiu  that  srere  then 
manipulated  to  derive  an  acoustic  subility  equation  in 
terms  of  the  more  readily  measurable  time  averaged  and 
fluctuating  quantities. 

The  developed  experimental  setup  dut  diis  dieoretical 
investigation  aims  to  describe  is  shown  in  Fig.  (1).  It 
consistt  of  a  15  meter  long,  3.75  cm  x  7.5  cm  duct  with 
the  multi-diffusion  flame  burner  (MDFB)  shown  in  Fig. 
(2)  installed  on  the  bottom  wall.  Though  initial 
experiments  will  be  performed  nsing  cold  flow,  the 
MDFB  can  be  used  to  simulate  the  behavior  of  a  gas 
phase  solid  propellant  flame.  The  MDFB  consisu  of 
1537  hypodermic  tubes  ananged  in  a  29  x  53  matrix. 
During  tet  flow  experiments,  oxidixer  flow  is  supplied 
through  the  hypodermic  tubes,  simulating  the  flow  of 
combustion  pr^ucts  from  burning  oxidizer  particles 
(e.g.,  ammonium  perchlorate)  in  a  composite  solid 
propellant,  and  the  fuel  is  supplied  through  the  spaces 
between  the  hypodermic  tubes,  simulating  die  flow  of 


pyrolysis  producu  from  the  fuel  binder,  h  cold  flow 
experiments,  room  temperature  air  is  supplied  through 
both  the  tubes  and  the  surrounding  spaces.  Tiro  acoustic 
drivers  attached  to  opposing  duct  walls  just  upstream  of 
the  exit  plane  are  us^  to  excite  a  stanttog  Imgitudinal 
acoustic  wave  in  the  duct,  simulating  an  axial  insubility 
fai  a  rocket  motor.  Row  is  injected  at  die  upstream  end  of 
die  duet  to  simulate  die  axial  core  flow  in  a  rocket  motor. 
This  head  end  injector  is  made  of  a  porous  plate  that 
behaves  as  an  acoustically  Tiard'  termination.  The 
location  of  the  MDFB  relative  to  the  acoustic  field  can  be 
changed  by  axial  translation  of  the  porous  plate. 
Directly  above  the  MDFB.  die  the  flow  through  the 
burner  turns  in  the  direction  of  the  core  flow  (i.e..  along 
the  axis  of  the  duct). .  This  is  where  flow  turning  losses 
are  expected  to  be  observed.  (>ptical  windows  m  this 
region  allow  LDV  measuremenu  of  acoustic  and  mean 
velocities  and  a  water  cooled  pressure  probe  allows 
acoustic  pressure  measurement. 

To  determine  the  behavior  and  relative  orders  of 
magnitude  of  the  various  urms  that  appear  in  the 
developed  subility  expression,  a  simple  numerical 
simulation  of  the  experimental  setup  that  will  be  used  to 
study  the  flow  turning  loss  was  carried  out.  The 
simulation  required  the  derivation  and  numerical  solution 
of  the  perturbed  mass  and  axial  momentum  conservation 
equations  for  the  proposed  experimental  setup.  The 
derived  solutions  were  substituted  into  the  developed 
acoustic  subility  equation  to  determine  the  relative 
orders  of  magnitude  of  the  various  Urms  that  sppear  in 
die  equation  and  whether  the  numerical  solution  satisfies 
the  developed  expression  for  the  growth  rau  o  of  the 
oscillations  in  the  "forced*  experiment.  These 
calculations  were  carried  out  for  a  number  of  different 
cxperimenul  configurations  and  operating  conditions, 
and  in  all  cates  excellent  agreement  between  the  results 
of  the  numerical  simulation  and  the  developed  subility 
equation  was  obtained. 
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MULTI-ISFFUSIC»4  FLAME  BURNER 
Rgure  A  schematic  of  the  experimental  setup  used  in  the  flow  turning  mvestigation. 
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The  initia]  goal  of  dda  iBvaad|atioa  «u  lo 
analydeall^  varify  dia  wfcianca  of  iha  flow  atmiag  loaa 
fey  dariviiig  an  oqnation  for  iha  aeoastic  tUbility  of  • 
ooBibiuior  fey  aa  approach  difTcmt  than  diat  aaod  fey 
Ctdick.  fli  order  to  make  the  iavaatigation  as  sdaight 
forward  at  possible,  and  make  it  qiplicable  to  the 
proposed  axperinenul  setup  (see  Fig-  1),  the 
favesitgaiad  geonstiy  coosisied  ^  a  2>D  raeiangular  duet 
afhsi^Handlcai^L.  The  fbmuilatien  is  based  upon 

a  procodure  developed  by  Canoell  and  Hai(^.  It  begias 
widi  oonsiderations  of  die  mass,  momeanan,  and  energy 
eonservation  aviations.  The  ntOiaad  two  diniaasioaal 
■MSS  and  nuimennun  eensarvadon  aquations  an: 


Rgun  2.  Schematic  of  the  inn1ti<diRi»ion  flame  burner 
■sed  to  simulau  combustion  of  composite  solid 
propellant  Aids  during  flow  ouiiiag  experimentt. 


whan  p  is  dM  density,  u  and  v  arc  the  axial  and  vardcal 
mmpmatu  of  the  velociqf  vector  V,  p  is  the  pressure 
and  m  is  dM  mass  flux  vector  pV.  Ihis  perdeular  form  of 
the  momentum  aquation  has  been  chosen  for  iu 
convcaieaoc  later  fa  the  analysb.  Assumiitg  An  ^ 
flew  is  iseanopic.  Eq.  (2)  bacomes 


udMn  hf  is  dM  total  anihafay  per  unit  mast.  At  a  result 
of  the  imnopic  flow  assumption,  this  analysiB  cannot 
account  for  losses  due  to  particle  damping  and  heat 
traiufer,  or  gas  phase  heat  addition  combustion 
processes,  h  a  more  general  treatment,  an  attempt  must 
be  made  to  consider  aon-iMeatropie  aflecu  due,  for 
example,  to  heat  mnsfer  and  viscous  processes. 

The  developincnt  of  the  acousdc  subiliiy  equation 
by  an  energy  balance  qiproach  staru  with  the  foUowfag 
form  of  the  energy  equation 

fli'  (4) 

Sfacc  die  acoustic  energy  is  a  second  order  acoustic 
quantity,  dl  first  and  second  order  terms  must  be  retained 
fa  the  andysis  while  third  order  and  higher  order  teems 
may  be  neglected.  Thus,  Eq.  4  is  expanded  to  second 
otde  and  time  averaged  to  yield 

•(P-  (mshTo)+V-  (mi  hn)*^-  (mofan)]k  (g) 

where  dM  quantities  inside  a  feradcct  o  are  averaged  over 
dM  subscripted  variable  and  the  numcried  subscripu 
fadieaie  the  order  of  term.  To  eliminate  the  dependence 
upon  y,  Eq.  (S)  is  averaged  over  the  height  (iM..  dw  y 
directiw)  of  constam  area  duet  to  obtafa 

^(pbT-pltJ)^^  -  -^(((P«)i  hn))i.. 


*^(((P  »>  M+dP  »>  l»n)li.> 


♦i{[UP»)^^o)+ftP»X**n)+((P'')bWj)S)i 

whan  H  is  tfie  height  of  Ae  ducL 

EqoatioB  (6)  includes  second  order  terms  whose 
sohiiioo  ere  generally  not  known  end  difficult  to  obtain. 
It  is  possible,  however,  to  use  Eqs.  (1)  end  (3)  to 
aliininate  dwse  tenns  in  Eq.  (6).  In  order  to  accomplish 
this,  some  manipulationi  of  ^s.  (1)  and  (3)  must  be 
performed.  WiA  this  in  mind,  die  seroA  order  part  of 
eq.  0)  (ije..  die  steady  i^ete  energy  equadon)  is  obtained 
by  retaining  an  terms  that  are  independent  of  time  and 
neglecting  terms  of  second  order  and  lugher  in  die  mean 
Mach  mtmber,  resulting  in  the  foUowmg  equation: 

Vhj5  -  0  (7) 

Neat,  Ae  continuity  equation,  Eq.  (I).  was 
multiplied  by  Ah<fQ  and  the  resulting  equation  added  to 
die  product  of  Eq.  0)  and  m.  The  resulting  equation  was 
then  averaged  over  y  and  t  to  obtain  Ae  following 
equation  where  all  terms  of  the  order  of  Mach  number 
aquwe  or  higher  ware  ruglected: 

^(((P«)i>»lo)|v,  ♦  ^(C((P 


AnoAer  needed  equation  is  obtained  by  forming  the 
dot  product  of  the  second  order  momentum  aquation,  Eq. 
(3),  wiA  Ae  aeroA  order  matt  flux  to  obtain  Ae 
following  expression  where  all  terrru  of  Ae  order  of 
Mach  number  square  or  higher  were  neglected: 

Equation  (9)  is  now  everaged  aerou  Ac  duet  to  yield 


(10) 


Equations  (S)  and  (10).  while  not  physically 
meaningAil  A  diem^as.  ate  fe  forms  that  are  useful  for 
eimpliiying  Ae  energy  equation.  Eq.  (6).  Substitution 
Eqs.  (8)  end  (10)  Ato  Eq.  (6),  rearranging  and 
aimplify^  the  resulting  etqirettion,  and  integration  of 
die  result  over  Ae  lengA  of  the  duct,  Ae  following 
equatioo  A  obtaAed: 


Po«^^PlVo«l 

2  J 


where  L  and  H  are  Ae  lengA  and  height  of  die  duct  and 

A  the  xcroA  order  speed  of  sound.  Equation  (11)  A  Ae 
Attgral  form  of  e  conservation  equation  for  the  quantity 
A  Ae  time  derivative.  It  must  U  made  clear  Aat  AA 
quantity.  wAch  hu  uniu  of  energy,  does  not  rqiresent 
the  totri  acoustic  energy  of  the  system.  The  lack  of  a 
simple  physical  interpretation  of  Ait  quantity  A  Ae 
necessary  result  of  the  cancellation  that  resulted  from  the 
above  described  simplifications  of  Eq.  (6). 

LettAg,  for  convenience,  the  energy  Ategral  A  Ae 
above  Ategral  equal  E^;  that  A, 


(jL+esi+ELHifilda 

Upoi?  2  ' 


(12) 
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and  noting  Aat  the  growA  rate  of  each  of  the  acoustic 
energy  terms  A  Ae  tame  at  the  growA  rate  of  E^,  then 
the  exponential  growA  rau  o  can  be  obtained  from  Ae 


fBOowiac  aipwtrioB; 


2a  ■ 


at  ax 


(17) 


(13) 


^-0  (X-O.L) 
dx 


(It) 


«*«•  V  •  MBond  odv  aeouitk  qmaiiqr.  UxfagBqk 
(12)  ad  (13X  die  foOowjiig  •eonetic  stibUi^  equeiioB 
ii  afcgined  fteia  Eq.  (11): 


2aE*  -  .^|pi  •, 

,t\‘\ 


■  J^({“1 1|  C  «l  if  )«lx  ) 


(14) 


when  b  die  meei  flux  (pv) ««  of  die  weUe  (ediere  ell 
burning  b  eesumed  to  teke  pleaX  ediicli  ie  by  definiiim 
^itive  in  the  poeitive  y  direction.  Tbit  equation  ii  by 
biiani  (and  not  by  my  t^ioua  choice  of  sroiqnnga)  of  e 
form  tto  allows  comperison  with  the  acoustic  sti^qr 
aquation  developed  by  Culick^.  Culick's  oquatkm. 
simplifled  by  removing  the  terms  that  contain  the  effecn 
of  particulate  matter,  gas  phase  combustion,  and  cross* 
duct  temperature  gradious,  b 


(15) 


of  Eqi. 


(14)  and  (IS)  teveab  diat  the 


appears  only  in  Bq.  (14).  Thb  may  be  attributed  to  die 
^^lication  of  the  boundary  condidoiu  shown  in  Eq. 
(It)  in  the  solution  of  Eq.  (16).  The  boundary 
conditions- force  the  approximate  vadocity  ampUtude  to 
aero  at  x  ■  0  and  x  ■  L.  In  thb  case,  die  omit^  term  b 
negligible.  The  benefit  of  the  acoustic  stability  equation 
givm  in  Eq.  (14)  b  that  it  b  appUcable  to  more  general 
types  of  oscillations.  The  term  associated  widi  the  flow 
ttining  loes  b  also  praecnt  in  Eq.  (14);  dial  b. 


nOWTUlNlNOLOSSn 


(19) 


hspection  of  thb  terra  thows  that  it  b  negative  if  mast 
mters  the  system  dirough  die  top  or  bottom  boundary, 
implying  diat  it  tends  to  reduce  a  and,  dius,  suMbe  the 
aystam.  On  die  other  hand,  if  maw  exiu  the  qrstem 
t^ugh  the  top  or  bottom  boundary  the  term  b  idways 
positive,  udiich  increases  a  md  destebilites  the  system. 
Unfortunately,  a  one*dimentionaI  analysis  caiuiot 
provide  buight  teto  the  detaib  of  the  flow  turning  lost 
mechanismt. 


Ilib  analysb  has  confirmed  diat  the  flow  turning 
low  b  a  part  of  the  one  diinensional  acoustic  subility 
aquation.  A  secondary  result  of  thb  malysb  ww  a 
confirmation  of  Vm  Moorhera’s^  conclusion  that  the 
flow  taming  tarm  b  not  absent  flora  a  three  dimenskmal 

malysb  as  proposed  by  Culick^,  but  that  the  term  cm 
only  be  isolated  w  a  sbigb  term  upon  averaging  over  the 

crow  section  die  duct  As  suggwted  by  Baum'*  die  flow 
turning  term  describes  a  combination  of  several 
processes  diet  result  whm  flow  b  injected  transverse  to 
the  acoustic  motions.  The  habBiqf  to  isolate  a  flow 
toning  term  fat  a  three-dimensional  analysis  suggesu 
dial  dmeloping  m  viderstmding  of  the  detaib  of  the 
mechanism  that  centrob  die  flow  toning  low  mqr  be 
vary  dUBcab. 


where  *  taprasentt  the  ampbude  of  a  harmonic  flmetim 
and  the  sobwript  I  signifiw  approximate  aohitions  for 
foe  acpttsric  pressure  and  velod^,  udiidi  are  solutioiis  of 
foe  xlusicd  acoustic  problem  governed  by  die 
following  aquadons  and  boundary  conditions: 


(16) 


BanimbtiiMi  far  die  EMerinuailel  Smdv 

The  one  dimensional  eubOify  aquation  shown  u 
Eq.  (14)  ww  derived  by  expanding  the  solutions  in  terms 
of  first  md  second  order  qumtities.  While  such  a 
procedure  b  theoretically  sound,  b  b  hard  to  oorrelaie 
Srst  md  order  qumtitiw  widi  rnOwured  experimmul 
data.  Therefore,  a  new  theoretical  formidation  that  cm 
be  used  in  the  experimmtal  determination  of  the  flow 
toning  low  hw  bem  developed  and  it  b  presented  in 
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dUi  Mctian.  Starting  with  tht  fbltowing  forms  of  dw 
■was  cantinuity  and  nial  momentum  acpiations 


(20) 

p»L+p«»L+p»il+§E.O 
*^at  ax  ay  ax 

(21) 

and  o^^rcssing  oadi  dependent  variable  as  a  sum  of  a 
tima  averaged  eonqmnent  and  a  fluctuating  component 
one  obtains 


•♦■~(pV+p  v*+p*V+p*  v’)«  0 

iy  (22) 


(23) 


when  the  bar  and  a  prime  denote  a  time  averaged  and  a 
fluctuating  quantity,  nspeciively.  Neglecting  producte 
of  perturbations  (ndiich  nil!  be  slwwn  to  be  equivalent  to 
Mglacting  temu  of  third  order  or  higher  in  the  result) 
and  subtracting  the  time-invariant  forms  of  the 
equations,  Eqs.  02)  and  (23)  become 


^♦^(P«’+P*»)+^(P»*+P'3)-  0 
dt  dx  dy  (24) 

Py^+P  ^(Brf)+P*3|^+pv|il 
ac  9X  ax  ay  ^ 


♦p  V*— ♦p'v— 0 
dy  dy  dx 


(25) 


Aasuming  that  the  fletw  b  bemropic.  as  was  done  in  dw 
ptnvious  analysb,  mult^lymg  ^s.  (24)  and  (25)  by 


_  ,  nspcctively,  adding  the 

nsuhiiig  equadflcis,  and  averaging  the  resulting  equation 
;  die  duct  (i-On  y)  oew  obtains: 


i- 


P^  *  *  I  f 


Using  die  relationship  given  in  Eq.  (13X  Eq-  (26)  can  be 
used  to  obtain  die  following  expreuion  for  the  complex 
growth  rate  constant  of  the  oscillations. 


Cotnparbon  of  Eqs.  (14)  and  (27)  reveab  that  the 
fem  of  dw  equations  b  idetuical.  However,  it  must  be 
noted  diat  die  quantities  involved  in  dwse  equations  are 
defined  differendy,  and  therefore  the  equations  are  not 
the  same.  Because  of  the  way  the  terms  heve  been 
defined,  Eq.  (27)  lends  itself  more  readily  to  the 
experimental  verification  of  die  flow  turning  losses. 

Numerical  Simulation 

In  order  to  determine  the  behavior  and  relative 
esders  of  magnitude  of  dw  various  terms  that  qipear  in 
the  developed  stability  expression,  Eq.  (27),  a  numerical 
simulation  of  the  experimental  setup  that  will  be  used  to 
study  the  flow  turning  loss  (see  Fig.  (1))  wu  carried  out 
The  simulation  involved  derivation  and  numerical 
soludm  of  the  perturbed  mass  and  axial  momentum 
conservation  aquations  for  the  proposed  experimanul 
aan^;  that  is. 
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Eqnatioiu  (29)  md  (30)  me  derived  etcwning  ihet 
1)  the  flow  ie  imiropic,  2)  the  pcrturbatioiu  are 
hinnoiiic  in  time.  3)  die  y>dqiaidenoe  of  the  dqicndeat 
verieUes  it  faiown,  and  4)  die  pretture  at  the  apetream 
hyodw  plme  it  known  (ia,  muturad). 

The  derived  matt  and  momentnm  oontervetion 
equationi  were  then  inugrated  along  die  exit  of  the 
experimental  actup  and  their  aolutiont  eubttituted  into 
the  developed  acouttie  tubility  equation,  Eq.  (27),  to 
determine  the  relative  ordert  of  maj^aide  of  die  varioiit 
termt  that  appear  in  the  aquation  and  whether  the 
numerical  tolutioru  tatitfy  Eq.  (27).  Since  the  tohidont 
tifflulaie  the  behavior  of  a  "forced"  experiment,  the 
amplitude  of  the  otcillaiiont  don't  fawreate  with  time  and 
the  growth  rate  a  must  be  tcro.  These  calculaiiont  were 
earried  out  for  a  number  of  different  experimental 
conflgurations  and  operating  eonditiont,  and  ta  all  eaaet 
excellent  agreement  between  die  retulu  of  the  numerical 
aifflulation  and  the  developed  tubiliqr  aquation  wet 
obtained. 

Retulu  of  a  typical  simulation  arc  shown  in  Figt. 
(3)  and  (4).  In  order  to  simplify  the  calculatioiu,  die 
attaittaneet  of  the  upstream  injector  wall  and  burner 
curface  were  utumed  to  be  to  xcro.  Furthermore,  it  was 
assumed  that  the  amplitude  of  the  pressure  oacQladon  at 
the  injector  face  is  200  Pa,  the  mean  pretture  is 
atmot^eric,  die  burner  it  10  cm  long  and  is  located  70 
cm  downstream  of  the  head  end,  and  that  the  mean  Madi 
number  of  the  ir^jected  flow  at  the  burner  surface  it 
0.004,  which  simulates  the  Mach  number  of  tha  flow  at 
the  multi-diffusion  flame  burner  surface  fai  a  cold  flow 
experiment.  Figure  (3)  describes  the  axial  variations  of 
the  amplitude  ^  pharo  of  the  acoustic  pressure  along 
die  experimental  setup.  The  location  of  the  burner  is 
shown  on  the  axis.  It  it  seen  that  for  the  investigated 
Mach  number,  the  effect  of  matt  addition  at  die  bmer 
surface  upon  the  pretture  it  small.  The  axial  variatiofu 
of  the  magnitudes  of  the  various  terms  that  appear  in  Eq. 
(27)  arc  described  in  Fig.  (4).  The  sum  of  these  tatins, 
whi^  equals  a.  it  alto  presented  and  it  equals  aero  at  aD 
locstiont,  in  agreement  with  the  theory.  It  thotdd  be 
noted,  however,  that  the  calculated  magnitude  of  die  flow 
turning  lost  for  the  experimental  setup  reveals  that  it  it 
quitt  small.  This  it  a  consequence  of  the  low  fayeetion 
Mach  number  and  the  small  perturbation  amplitude 
(compared  to  those  experienced  in  unsuble  solid  rocket 
motors),  and  does  not  suggest  that  the  flow  turning  lots 
k  always  small.  Rgure  (4)  alto  shows  diat  aD  the  tatms 


that  are  first  order  In  Mach  number  are  of  similar 
nugnitude.  Because  diese  tamu  arc  of  similar  order, 
Chen's^  experimental  approach,  whidi  did  not  account 
for  aU  relevant  terms  in  the  staUlity  aqpatioii,  xould  not 
affaedvefy  measure  the  flow  ttnning  kas. 

Tim  numerical  shnulation  eras  nsad  to  study  the 
affaett  of  several  parameters  upon  the  magnitude  of  the 
flow  turning  lost.  One  of  ftete  parameters  was  the 
location  of  the  burner  relative  to  the  standing  acoustic 
erave.  The  resulu  of  this  soidy  are  presented  in  Fig.  (S). 
It  indicates  that  the  flow  turning  lou  depends  upon 
where  the  injection  occurs  relative  to  t^  standing 
acoustic  wave.  The  (low  turning  loss  is  maxunum  and 
minimum  when  the  flow  is  injected  near  an  acoustic 
velocity  maximum  and  minimum,  respecdvely.  The 
second  parameter  studied  was  the  mean  core  flow  Mach 
number,  ediich  was  varied  from  M  ■  0.0  to  M  ■  0.1.  This 
series  of  simulations  showed  that  die  magnitude  of  the 
flow  turning  lots  was  independent  of  the  core  velocity 
over  the  invettigated  range.  Third,  it  was  demonstrated 
that  the  flow  turning  lots  varies  linearly  srith  the 
injection  Mach  number  over  an  invettigatad  range  of  M 
■  0.0  to  M  ■  0.1.  Finally,  the  dependence  of  the  flow 
Mtning  lou  i^on  the  amj^nide  of  the  acoustic  prcuure 
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Figure  3.  NumericaDy  predicted  amplinide  and  phase  of 
die  pteuute  otcillstion  k  a  duct  with  side  ityeetion  near  a 
pressure  minimum.  The  phau  is  meuuied  srith  respect 
to  the  pressure  signal  at  x  ■  0.  The  injection  occurs 
between  x  ■  0.7  and  x  ■  0.t. 
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OMiHalion  wm  iiivestigitsd  and  it  wu  fboad  that  it  it 
fropoitioiul  to  the  square  of  the  mplinide,  m  tacpected. 
Note,  however,  that  the  growth  rate  a  it  proportional  to 
die  Im  of  energy  divid^  by  the  total  toiergy,  and  to  a 
it  adcpcndent  of  the  amplitude  of  die  otefllationt. 

Concludine  Remetka 

The  inclusion  of  flow  turning  lottet  in  oodea  diet 
predict  the  sttbiliiy  of  solid  rocket  motors  hat  been  die 
subject  of  some  oontroverty.  The  theoretical  aialyait 
performed  in  diit  pqwr  hat  verified  the  axittenee  of  die 
term  and,  therefore,  aidi  in  validating  its  inclution  fa 
one  •dimensional  atability  predictions  modelt. 

Furthermore,  this  study  eonflrmed  Van  Moorhem's^ 
conchaion  that  the  one  dimensional  analysis  contains 
no  information  that  it  not  pretent  in  properly 


S  -THE  FLOW  TURNING  USS 

Rguic  4.  Numerical  prediction  of  the  relevant  terras  on  the 
fi^t  hand  side  of  the  1  •dimensional  acoustic  lubiliiy  equation 
for  a  duct  with  tide  injection  from  .7  to  J  meters  from  die  head 
and,  in  this  cate  near  a  pressure  minimum. 


performed,  three  dimensional  atability  analytes.  The 
flow  turning  term  cannot,  however,  be  isolated  as  a 
single  term  in  the  three  dimensional  analysis, 
suggesting  that  the  term  it  due  ib  an  artificial  grouping 
of  terms,  and  not  a  tingle  mechaiutm.  This  result  dwly 
saublithet  that  the  mechanitmt  favoNed  fa  the  'flow 
turning'  term  arc  already  accounted  for  in  a  multi¬ 
dimensional  analysis,  and  the  addition  of  the  term  to  a 
multi^dimensional  analysis  it  not  justified.  For  this 
reason,  die  study  of  fae  flow  turning  lots  as  it  is 
currently  understood  must  be  tettricted  to  a  one- 
dimensional  analysis,  r-. 

The  second  result  of  diis  theoretical  analysis  was 
the  development  of  an  acoustic  stability  aquation  that 
can  be  used  to  experimentally  determine  the  flow  tuming 
lots.  The  ttandsrd  theoretical  ^iproach  in  which 
solutions  arc  expanded  in  power  series  of  the  amplitude 
and  Mach  number  does  not  produce  resulu  that  can  be 
readily  used  in  the  interpreution  of  experimental  dau 
because  of  the  vague  eonnecUon  between  the  theoretical 
and  measured  variables.  To  resolve  this  problem,  a 
different  dieoretical  approach  in  which  the  dependent 
variables  are  written  at  the  sum  of  a  time^averaged  and  a 
fluctuating  component,  which  can  be  measured,  was 
developed.  This  model  together  with  a  simple  numerical 
simulation  of  the  proposed  experimental  set-up  were 
then  used  to  study  the  expected  behavior  of  the 
experimental  setup  and  relative  orders  of  magnitude  of 
the  various  terms  in  the  derived  subility  aquation.  This 
analysis  was  canied  out  for  a  number  of  different 
cxpetimenul  configurations  and  operating  conditioiu, 
and  in  all  cases  excellent  agreement  between  the  results 
of  the  numerical  simulation  and  the  developed  stability 
oquatioR  was  obtained. 

The  next  in  this  study  is  the  expcrimcnul 
investigation  of  the  flow  tuming  lost  using  the  newly 
developed  subility  equation.  Measurements  will  be 
performed  to  determine  the  validity  of  the  one^ 
dimensional  subility  equation.  The  study  will  include 
experimental  deurminsiion  of  the  cffecu  of  various 
parameters,  such  as  the  location  of  the  flow  turning' 
region  relative  to  the  standing  acoustic  wave,  the 
frequency  and  amplitude  of  the  osdllation.  and  the 
magnitu^  of  the  core  flow  and  the  injection  velocities, 
upon  the  magnitude  of  the  flow  tuming  loss.  These 
resulu  will  be  then  compared  with  analytical  and 
aumarical  rasnlu. 
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ABSTRACT 

This  paper  describes  an  experimental 
investigation  of  the  flow  turning  loss  in  unaahle  solid 
propellant  rocket  motors.  It  is  caused  by  the  interaction 
of  the  radially  directed  flow  of  combustion  products 
leaving  the  propellant  surflue  with  the  axial  acoustic 
field.  The  erqrerimental  setup  consisted  of  a  rectangular 
duct  with  side  injection  that  simulated  the  combustion 
products  leaving  the  propellant  surfice.  Acoustic 
driven  were  used  to  excite  oscillations  that  simulated 
an  axial  instability.  The  mean  and  acoustic  velocities 
and  the  acoustic  pressure  oscillations  were  measured  in 
a  region  where  the  flow  turning  loss  was  expected  to 
occur.  The  measured  data  were  used  to  evaluate  the 
terms  that  appear  in  an  acoustic  stability  equation 
containing  the  various  loss  and  gain  terms  iricluding  the 
flow  turning  loss.  Experiments  were  performed  to 
determine  the  effects  of  the  iiyection  and  mean  flow 
Mach  numbers,  the  amplitude  of  the  oscillation  and  the 
location  of  the  burner  relative  to  the  standing  wave 
upon  the  magnitude  of  the  flow  turning  loss  in  the 
investigated  control  volume.  The  flow  turning  loss  was 
shown  to  depend  strongly  upon  the  injection  velocity 
and  the  location  of  the  injection  relative  to  the  standing 
wave.  Measured  data  compared  well  with  numerical 
predictions. 

INTRODUCTION 

The  objective  of  this  investigation  is  to  develop  a 
better  understanding  of  the  flow  turning  loss  and  other 
processes  that  contribute  to  the  driving  and  damping  of 
axial  instabilities  in  solid  propellant  todeet  motors. 
This  nnderstaixiing  is  necessary  for  the  development  of 
piactical  methods  for  the  design  of  stable  solid 
propellant  rocket  motors. 
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The  occurrence  of  combustion  instability  is 
determined  Ity  the  relative  magnitudes  of  the  various 
driving  and  damping  mechanisms  within  the  rocket 
motor  that  add  and  remove  energy  from  the  oscillations, 
respectively.  ^  the  acoustic  energy  gains  outweigh  the 
loses,  smidl  perturbations  amplity  to  levels  that  can 
adversely  affect  guidance,  thrust  history,  and  the 
structural  integrity  of  the  motor'**.  In  order  to  eliminate 
or  reduce  the  severity  of  acoustic  oscillations  resulting 
from  combustion  instabilities  in  solid  rodeet  motors,  the 
ma^hanitmit  by  iduch  energy  is  added  to  and  removed 
from  these  oscillations  must  be  identified  and 
understood.  It  is  generally  acoqrted  that  the  interaction 
between  the  solid  propellent  combustion  process  and  the 
acoustic  oscillations  provides  the  energy  necessary  to 
initiate  and  maintain  instabilities  in  rodeet  motors. 
Other  fretors,  such  as  nozzle  damping,  viscous 
dissipation,  and  flow  turning  generally  remove  energy 
from  the  oscillations  and  thus  stabilize  the  system. 

The  velodty  of  combustion  products  is  generally 
radial  at  the  burning  solid  propellwt  surface.  As  these 
gasses  move  away  from  the  prt^Uant  and  are  entrained 
in  the  core  flow,  tbQ'  are  turned  from  radial  to  axial,  as 
shown  in  Fig.  1.  During  this  turning  process,  the 
combustion  products  that  leave  the  burning  propellant 
surftce  with  no  axial  acoustic  momentum  roust  acquire 
axial  acoustic  energy  firom  the  core  flow  osdUations. 
This  phenomena  was  first  discussed  by  Culick’  who 
argued  that  since  the  combustion  products  entering  the 
core  flow  acquire  acoustic  energy  from  the  axial 
acoustic  field  (i.e.,  the  core  flow  oscillations),  the  axial 
acoustic  fidd  loses  acoustic  energy,  which  tends  to 
stabilize  the  tystem. 

The  flow  turning  phenomenon  has  been  studied 
ejqrerimentally  by  Herseb*'^  and  Chen*  and  analytically 
by  Baum*-'*  and  Hegde"'”.  In  these  investigations,  the 
behavior  of  acoustic  waves  in  a  diannd  in  which  the 
adid  propellant  combustion  process  was  simulated  b>' 
mass  addition  through  a  si^  wall  was  investigated. 
Hersch’s  conclusions  were  inconsistent  with  the 
behavior  of  the  flow  turning  loss  as  predicted  by  Culick. 
Also,  the  accuracy  of  the  results  of  both  Hersch's  and 
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Figure  1.  Schematic  of  a  solid  propellant  rodcet  motor  showing  the  flow  of  combustion 
products  turning  from  radial  to  axial  as  they  are  entrained  into  the  core  flow  and  the 
pressure  distribution  characteristic  of  axial  acoustic  instability. 


Chen's  investigations  are  open  to  question  due  to  the 
methods  employed.  Baum'^-'*  computational  studies 
show  that  the  flowfield  near  the  wall  is  extremely 
complex,  involving  a  multitude  of  vortical  souctures 
and  energy  exchange  between  the  mean  flow  and  the 
oscillations.  Hegde  et  al"'”  used  a  linear  analysis  and 
hot  wire  measurements  to  determine  the  acoustic 
diaracteristics  and  losses  of  a  duct  containing  a 
standing  axial  acoustic  wave  and  mass  addition  through 
the  side  wall.  His  results,  in  contrast  with  those  of 
Baum,  suggest  that  a  linear  model  may  be  adequate  for 
describing  the  important  features  of  the  flow  fidd  in  an 
unstable  solid  propellant  motor. 

The  goal  of  this  project  is  to  determine  the  validity 
of  the  flow  turning  loss  term  that  is  included  in 
currently  used  one  dimensional  linear  stability 
prediction  codes.  The  first  step  was  to  verify  the 
existence  and  mathematical  form  of  the  flow  turning 
loss  as  defined  by  Culidc*.  To  this  end,  a 
ooe'dimensional  acoustic  stability  equation  was  derived 
using  an  energy  balance  approach  similar  to  that  used 
by  CantrdI  and  Hart'^  with  the  addition  of  flow  rotation 
effects.  This  anafysis  confirmed  that  the  classical  flow 
turning  loss  term  is  indeed  a  part  of  the  one 
dimensional  acoustic  stability  equation.  This  analysis 
also  cmfinned  Van  Moorhm'^’  conclusion  that  the 
flow  turning  loss  is  not  absent  from  three  dimensional 
suibility  models  as  suggested  by  Culidc'*.  It  was  drawn 
that  the  "classical*  flow  turning  loss  term  results  from 
the  cross  sectional  averaging  used  to  reduce  the  three 
diinensional  problem  to  one  dimension.  The  inability 


to  isolate  the  flow  turning  term  in  a  three  dimensional 
stability  equation  suggests  that  developing  an 
understanding  of  the  details  of  the  mechanisms 
involved  in  the  term  may  be  impossible. 

The  second  phase  of  this  program  focused  on  the 
develc^ment  of  an  acoustic  stability  equation  that  was 
better  suited  to  the  interpretation  of  erqreiimental  data. 
The  one  dimensional  acoustic  stability  equation  is 
ordinarily  derived  Iry  expanding  the  governing 
equations  in  terms  of  first  arid  second  order  quantities. 
W^Ie  such  a  procedure  is  theoretically  sound,  it  is 
difficult  to  oondate  first  and  second  order  quantities 
with  measured  data.  Therefore,  a  new  theoretical 
qiproacb  was  developed  which  simplifies  the 
interpretation  of  data  in  the  experimental  d^ermination 
of  the  flow  turning  loss.  This  approach  begins  from  the 
equations  of  mass  and  axial  momentum  conservation  in 
which  each  of  the  dependent  variables  is  erqiressed  in 
terms  of  a  time  averaged  and  a  fluctuating  component. 
Fm  example,  the  pressure  is  erqiressed  as 

F-T+f'  (11) 

where  ^  is  the  time  average  of  the  measured  pressure. 
The  procedure  is  eiqilained  more  thoroughly  in  Ref  17. 

The  equation  for  the  growth  rate  constant  deri\’ed 
by  this  method  for  a  2>D  rectangular  volume  of  length  I 
and  height /f  is 
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uteie  a  is  the  acoustic  growth  nte  constant,  £*  is  a 
term  with  units  of  energy,  v  and  v  represent  the  axial 
and  transverse  velocity  components,  respectively,  and 
1R4  is  the  mass  flux  at  the  wall.  The  terms  in  this 
equation  were  grouped  in  a  way  that  permits 
conqrarison  with  Culi^'s  result  As  noted  by  Culick’. 
the  proper  choice  of  groupings  is  not  obvious,  as  qxttial 
averaging  of  terms  makes  it  difBcuIt  to  physically 
interpret  these  terms.  The  term  in  Eq.  1  that  describes 
the  flow  turning  loss  is  given  by 


FLOW  TURNING  LOSS  - 

Inflection  of  this  term  shows  that  it  will  always  be 
negative  in  a  rodcet  motor  where  mass  enters  th^gh 


the  burning  surfaces.  For  such  a  configuration  this 
term  tends  to  reduce  a.  which  stabilizes  the  system.  On 
the  other  hand,  if  mass  exits  the  tystem  through  a  side 
wall  (e.g.,  as  through  a  T*bumer  nozzle),  this  term  is 
positive,  which  increases  a  and  destahilires  the  tystem. 

EXPERIMENTAL  STUDY 

The  experimental  study  was  conducted  in  the  setup 
shown  in  Fig.  2.  It  was  designed  to  simulate  flow 
phenomena  near  the  walls  of  a  solid  propellant  rocket 
experiencing  an  axial  acoustic  instability.  The  facility 
consists  of  a  duct  7.6  cm  (3.0  in.)  high  by  3.8  cm  (1.5 
in.)  wide  fay  3  m  (9.1  ft.)  long  steel  duct  with  a  flat 
flame  burner  mounted  on  the  bottom  wall.  The 
multi-difliision  flame  burner,  which  was  designed  to 
simulate  a  gas  phase  solid  propellant  flame,  consists  of 
a  matrix  of  hypodermic  tubes  which  supply  oxidizer 
flow,  simulating  the  flow  of  combustion  products  from 
burning  ammonium  perchlorate  particles.  Fuel  is 
supplied  through  the  q>aces  between  the  tubes,  which 
simulates  the  flow  of  pyrolysis  products  from  the  fuel 
binder.  In  cold  flow  experiments,  room  temperature  air 
is  injected  through  the  burner.  To  date,  only  cold  flow- 
experiments  have  been  performed.  Two  acoustic 
drivers  mounted  on  opposing  duct  walls  just  upstream 
of  the  exit  plane  are  used  to  excite  a  standing  acoustic 
wave  in  the  duct  that  simulates  an  axial  insubility  in  a 
rocket  motor.  Axial  core  flow  is  simulated  by  the 
injection  of  air  at  the  upstream  end  of  the  duct  through 
a  movable  porous  plate  injector.  The  location  of  the 
wall  mounted  burner  (and,  therefore,  the  flow  turning 


Figure  2.  Schematic  of  the  experimental  setup  used  in  the  flow  turning  studies. 
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legion)  relative  to  the  ttanding  acoustic  wave  can  be 
varied  by  axial  translation  of  the  core  flow  iiyector. 
Velocity  measurements  were  performed  using  an  LDV 
tytttm  through  quartz  windows  in  the  aide  walls 
adjacent  to  the  wall  mounted  burner.  Acoustic 
piessuies  were  measured  with  a  1/4  in.  BrOel  A  Kicr 
anicrophone  probe. 

In  a  typical  experiment,  the  components  of  the 
mean  and  acoustic  velocities  and  the  acoustic  pressure 
were  measured  on  a  rectangular  grid  in  the  region 
idrove  the  burner  where  flow  turning  loss  was  expected 
to  occur.  A  iypical  mean  velociQ'  field  is  shown  in  Fig. 


3.  In  this  configuration  the  core  flow  was  injected  from 
the  right,  the  burner  was  on  the  bottom,  and  the 
acoustic  drivers  were  located  downstream  (to  the  left)  of 
the  displayed  control  volume.  The  acoustic  intensity 
vector  field  for  this  configuration,  shown  in  Fig  4, 
reveals  that  for  these  test  conditions  (i.e.,  frequency, 
core  flow  velocity,  burner  iiyection  velocity,  and 
position  of  injection  point  with  respect  to  the  standing 
axial  acoustic  wave)  a  relatively  large  amount  of  the 
acoustic  energy  is  Erected  into  the  burner  where  it  is 
most  likely  absorbed.  This  turning  of  the  acoustic 
intensity  flow  represents  a  zeroth  order  (in  Mach 
number)  acoustic  energy  loss.  The  amount  ot  acoustic 


Figure  3.  Vector  plot  of  a  typical  measured  mean  velocity  field  in  the  flow  turning  region. 


Figured.  Diagram  of  typical  measured  acoustic  intensity.  <p'ir>. 
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iflieiisity  ■tombed  by  the  buiiier  dq)ends  190B  the  lest 
coodition. 

The  lesultt  of  two  ■cries  of  tests  are  presented  in 
F^.  S  and  6.  Figure  5  describes  the  results  of  a  study 
of  die  dqiendence  of  die  Ikiw  turning  loss  upon  the 
■ean  core  flow  Mach  number  h^.  These  data  were 
■totained  by  varying  the  mean  flow  rate  iiyected  through 
die  porous  iryector  at  the  entrance  of  the  mid 
Theoretical  analysis  and  numerical  predictions  indicate 
diat  die  flow  turning  loss  does  not  depend  upon 
which  contradicts  the  mqieriinental  results.  This 
discrepancy  results  from  the  &ct  that  the  numerical  and 
dworetical  studies  investigated  the  behavior  of  a  control 
volume  that  qianned  the  entire  cross  section  of  the  duct, 
flircing  all  the  injected  flow  to  turn  in  the  direction  of 
die  axis.  In  contrast,  the  experimental  study 
investigated  the  flow  behavior  within  a  control  volume 
that  did  not  extend  from  the  burner  tatfut  to  the  top 
wall.  Gmsequently,  as  shown  in  Fig.  3.  the  flow  at  the 
•op  of  the  control  volume  is  not  "fully  turned*  into  the 

direction.  As  the  core  flow  Mach  mimber 
ftmreases.  the  amount  by  which  the  flow  injected 
through  the  wall  mounted  burner  is  turned  also 
incteases.  Therefore,  it  is  expected  that  as  the  core  flow 
Mach  nuniber  increases,  the  measured  flow  turning  loss 
will  also  increase,  which  is  consistent  with  the  trend 
exhibited  by  the  data  presented  ia  Fig.  S. 
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Figure  6.  Axial  variation  of  the  flow  turning  loss  term 
for  3  burner  iryection  Mach  numbers. 

this  relationship  is  not  linear.  The  discrqiancy  is 
again  believed  to  be  due  to  the  incomplete  flow  turning 
in  the  investigated  control  volume.  As  M^  increases, 
the  flow  turning  region  diickens  and  the  flow  turning 
process  requires  a  larger  region  above  the  burner  for 
completion.  It  is  expected,  therefor,  that  the  rate  at 
which  the  flow  turning  loss  in  the  investigated  control 
volume  increases  with  will  be  less  than  linear,  and 
this  agrees  with  the  observed  resultt. 

The  term  Ont  describes  the  contribution  of  the 
flow  turning  loss  to  the  growth  rate  constant  a  of  the 
investipted  region.  The  eflects  of  the  investipted 
parameters  upon  Ont  are  shown  in  Fip.  7*13.  Figures 
7  and  S  show  the  effect  of  burner  iryection  velocity 
qx>n  Oni  in  the  investipted  region  for  two 
experiments.  The  results  shown  in  Fig.  7  were  obtained 
ftm  the  experimental  setup  configured  as  followa.  the 
flequenqy  of  oscillation  was  S50  Hz.,  the  mean  core 
flow  velocity  was  0.4  m/s,  and  the  irye^on  occurred  at 
a  pressure  node.  The  measured  Om.  varies  linearly 
with  the  mean  burner  iryection  velocity.  Figure  8 


Figures.  Axial  variation  of  the  flow  turning  loss  for 
3  core  flow  Mach  numbers. 

The  resultt  of  a  atu4y  of  the  dqiendence  of  the 
flow  turning  loss  apon  the  burner  irijection  Mach 
■umber  are  presented  in  Fig.  6.  In  this  case,  the 
flow  taming  loss  is  plotted  for  three  tests  in  which  tmly 
the  rale  of  mass  iry^on  ww  varied.  The  theoretical 
predictions  indicate  a  linear  relationship  between  the 
burner  irtjection  Madi  nuniber  and  the  magnitude  of  the 
flow  turning  loss.  While  the  experimental  results 
clearly  show  that  the  flow  turning  less  increases  arith 


Figure?.  Dependence  of  Onv^pon  the  mean  burner 
iiyeciion  Mach  number. 
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Aows  itwlts  ebtained  with  •  Bean  ooit  flow  vdociQr 
of  1.T7  n/$,  a  SSO  Hz.  oacillation,  and  the  flow  turning 
Rgion  located  between  a  piessnie  node  and  antinode. 
The  ateasuied  is  again  observed  to  vary  hoeariy 
with  the  burner  injection  velodQr.  This  linear 
sriatiofidiip  agrees  with  dw  theoretical  — MiDeriGal 
predictions. 


Figures.  Dependence  of  upon  the  mean  burner 
iqiection  Mach  number. 

Figures  9  and  10  show  the  dqiendenoe  of  ot^  in 
die  investigated  li^on  upon  the  mean  core  flow 
vdodty;  the  data  in  Fig.  9  were  measured  in  tests  with 
a  ftequenQ'  of  550  Hz.,  a  burner  iqjection  vdodty  of 
0.156  m/t,  and  the  burner  located  at  a  pressure  node; 
the  data  in  Fig.  10  were  dnairted  with  a  ftequency  of 
480  Hz.,  a  burner  irtjection  vdodty  of  1.02  m/S  and  the 
flow  tureing  region  located  between  a  pressure  node 
and  a  pressure  antinode.  In  both  aeries  of  tests,  the 
measured  o„t,  increased  with  dw  mean  core  flow 
vdodty  •rrd  appears  m  approadi  a  liiniL  Again,  as 
discussed  before,  it  is  believed  that  this  behavior  is  Are 
to  inoonqrlete  flow  turning  in  the  investigated  region. 
Numeried  simulations  predict  this  behavior  nAen  this 


KgurelO.  D^endenoe  of  ovrttqMn  the  mean  core 
flow  Mach  number. 

inoonqrlete  flow  turning  process  is  accounted  for.  It 
diould  be  pointed  out,  however,  that  the  difference  in 
die  trends  exhibited  by  the  data  in  Figs.  9  and  10  is  due 
to  the  fret  that  in  the  tests  described  in  Fig.  10  the 
magnitude  of  the  core  flow  vdodty  was  high  enough 
diat  its  effect  upon  the  flow  turning  loss  is  small. 

The  dqrendence  of  Om.  iq>on  the  amplitude  of  the 
excited  acoustic  wave  Is  presented  in  Figs.  11  and  12. 
Figure  11  describes  tests  conducted  with  mean  core 
flow  vdodty  of  2.25  m/s,  mean  burner  iryection 
vdodty  of  0.76  m/s,  300  Hz.  ftequency,  and  iryection 
between  an  acoustic  pressure  node  and  an  antinode 
while  Fig.  12  describes  teste  conducted  with  a  mean 
core  flow  vdodty  of  1.45  m/s,  mean  burner  injection 
vdodty  of  0.62  m/s,  a  550  Hz.  ftequenqr,  and  iigection 
at  an  acoustic  pressure  antinode.  The  pressure 
amplitudes  were  normalized  with  reqpect  to  the  largest 
anqrlitude  oscillation.  Both  flgures  diow  that  ttm.  is 
independent  of  the  aiq;>litude  of  osdllstion,  wtich  is 


Figure  9.  Dependence  of  0^  upon  fre  mean 
core  flow  Madmumber. 
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Figure  11.  Dyendenceofg^npontherron- 
dimensjond  aaqilitude  of  oscillation. 


Figure  12.  DqKodenoe  of  ttm  iqwn  the  non- 
dimensional  amplitude  of  oscillation. 

expected  in  a  linear  analysis  and  is  in  agreement  with 
tbraiy  and  calculations. 

The  dependence  of  the  measured  value  of  Oni 
qwn  the  location  of  the  flow  turning  region  relative  to 
tk  standing  acoustic  wave  is  shown  in  Fig.  13.  This 
cjqieriment  was  performed  with  a  mean  core  velodiy  of 
2.86  in/s.  mean  burner  iiyection  velocity  of  1.34  m/s. 
and  a  frequenor  of  550  Hz.  The  measured  Om.  ^ 
greatest  when  the  flow  turning  occurred  at  an  acoustic 
pressure  node,  and  the  smallest  when  the  flow  turning 
region  was  located  at  a  pressure  antinode.  This  tread  is 
in  agreement  with  the  theoretical  prediction  that  a^v  is 
proportional  to 


Rgurel3.  DqpendenoeofO(ni.tqx>n  the  location  of 
the  flow  ir^lection  point  relative  to  the 
standing  acoustic  wave. 

A  munerical  Emulation  of  the  investigated 
experimental  setup  using  measured  boundary  conditions 
was  performed  to  determine  whether  the  mqierimental 


results  could  be  modelled.  The  simulation  involved  the 
derivation  and  numerical  solution  of  the  perturbed  mass 
and  axial  momentum  conservation  equations  for  the 
developed  experimental  seuq>.  assuming  that  the  flow 
was  iaentropic  and  the  perturbatitms  were  harmonic  in 
time.  The  equations  were  Integrated  along  the  axis  of 
the  er^rimental  setup  and  their  solutions  were  used  to 
determine  the  values  of  the  various  terms  that  appeu  in 
the  acoustic  stability  equation.  The  resulu  of  a  typical 
simulation  are  shown  in  Figs.  14-16.  The  cross-duct 
average  of  the  measured  axial  acoustic  velocity  and  the 
ooneqxmding  calculated  values  for  a  typical 
experimental  configuration  are  shown  in  Fig.  14.  The 
agreement  between  the  measured  and  calculated 
velocities  is  excellent 


Figure  14.  Comparison  of  measured  and  calculated 
values  of  the  cross-sectional  average  of 
the  axial  acoustic  velocity  in  the  flow 
turning  region. 

Figure  15  presents  comparisons  between  the 
aNasured  and  calculated  acoustic  pressure  amplitudes 
and  phases  between  the  pressure  and  velocity 
oscillations.  Good  agreement  between  the  calculated 
pressure  and  measured  pressure  amplitudes  is  obsened, 
but  the  error  between  the  calculated  and  measured 
phase  values  is  large. 

The  calculated  and  measured  values  of  the  flov. 
taming  loss  for  this  test  configuration  are  shown  in  Fig. 
16,  which  shows  excellent  agreement  between  the 
aaeasured  and  calculated  values.  Note  that  the  value  of 
the  flow  trrming  loss  term  in  this  erqieriment  is 
negative,  indicating  a  gain  of  acoustic  energy.  This  is 
observed  because  of  the  nature  of  the  flow  in  the 
investigated  control  volume.  The  flow  iiyected  through 
the  burner  is  turned  into  the  axial  direction,  while 
simultaneously  the  core  flow  entering  the  control 
volume  from  upstream  is  being  turned  upward.  For 
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Fifure  15.  Comparisons  between  measured  and  calculated  cross^ectional  avenges  of  the  acoustic  pressure 
amplitudes  and  phase  angles  between  the  pressure  and  the  axial  velocity  oscillations  in  the  flow 
turning  region. 


- CALCULATED 

•  MEASURED 


Figure  16.  ConqMrison  of  measured  and  calculated  values  of  the  axial  variations  of  the  nte  of  flow  turning 
loss  across  the  investigated  centred  vtdume. 


aome  mqieiiniental  oonflguralions,  dw  flow  turning  loss 
due  to  the  turning  of  the  iiyected  flow  is  smaller  than 
dbe  Tflow  turning  gain”  due  to  the  turning  of  the  con 
flow.  This  flow  turning  gain  effect  is  Observed  only 
because  of  the  Aoioe  of  control  volumes;  h  will  not 
•oenr  in  a  volume  extending  ftom  waU  to  wall  wboe  a 
flow  taming  loss  will  alwi^  occur. 

CONCLUDING  REMARKS 

This  paper  iqxNis  the  icsults  of  an  experimental 
invastigatioa  of  Ite  flow  turning  loss  in  a  setup 


designed  to  simulate  a  regjoD  near  a  burning  propellant 
surfhoe  in  a  solid  propellam  motor  oqteriencing  an 
axial  instabili^.  The  mean  and  acoustic  velocities  as 
well  as  the  acoustic  pressures  were  measured 
throughout  a  oontrol  volume  encompassing  the  flow 
taming  region,  and  the  measured  data  were  used  to 
evaluate  the  flow  turning  loss  and  rdated  terms  in  an 
acoustic  stN>i%  equatkm  developed  eariier  under  this 
program.  These  measurements  diowed  that  the 
magnitude  of  the  flow  tumiag  loss  is  mainly  determined 
by  die  wall  ii|jection  velocity  and  die  lo^od  of  the 
injection  point  relative  to  the  standing  acoustic  wave. 
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This  ooocIusioD  agrees  wtth  die  pradictioDs  uedsting  6) 
dieoretical  models. 

Measured  boundafy  oonditions  were  used  in  a  7) 
pieviously  developed  mimerical  Emulation  of  the 
experimental  aetiq;>.  The  calculated  values  of  the 
cncs-aection  average  of  the  axial  acoustic  velocity.  8) 
acoustic  pressure  amplitude,  the  phase  between  the 
axial  velocity  and  pressure,  and  the  axial  variation  of 
the  flow  turning  term  were  compared  to  measured 
values.  The  inviscid  equations  uted  in  the  munerical  9) 
simulation  were  tf>le  to  adequately  predict  the 
amgnitudes  of  the  pressure  and  axial  velocity  variations 
(and,  therefore,  the  flow  turning  loss  term),  but  did  not  10) 
accurately  predict  the  phase  between  them.  Thisfiulure 
to  predia  the  zeroth  order  (in  Mach  number)  acoustic 
energy  loss  suggests  the  possibility  that  viscous  driving 
or  danq>ing  mechanisms  that  are  not  treated  in  this 
analysis  or  accounted  for  in  standard  stability  prediction  11) 
methods,  may  be  inqwrtant  In  order  to  investigate  this 
possibility,  detailed  numerical  modelling  of  the 
experimental  setup  using  a  full  Navier*Stokes  solver  is 
currently  underway.  Results  fiom  the  study  will  be  used 
in  conjunction  with  previous  results  to  determine  the 
inqxMtanoe  of  the  viscous  terms  in  axial  stability 
predictions. 
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